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fi 要 结构 方程 模型 是 心理 学 、 管 理学 、 社 会 学 等 学 科 中 重要 的 统计 工具 之 一 。 然 而 ,大 量 使 用 结构 方程 
模型 的 研究 忽视 了 对 该 方法 的 统计 检验 力 进行 必要 的 分 析 和 报告 ， 在 一 定 程 度 上 降低 了 这 些 研 究 的 结果 的 证 
明 效 力 。 结 构 方程 模型 的 统计 检验 力 分 析 方 法 主要 有 Satorra-Saris 法 、MacCallum 法 与 Monte Carlo 法 三 类 。 
其 中 Satorra-Saris 法 适用 于 备 择 模 型 清晰 、 检 验 对 象 相对 简单 、 检 验方 法 基于 分布 的 情形 ; MacCallum 法 适 
用 于 基于 Xe 分布 的 模型 拟 合 检验 且 备 择 模 型 不 明 的 情形 ; Monte Carlo 法 适用 于 检验 对 象 相 对 复杂 、 采 用 模拟 
或 重 抽 样 方法 进行 检验 的 情形 。 在 实际 应 用 中 ,研究 者 应 当 首先 判断 检验 的 目的 、 方 法 以 及 是 否 有 明确 的 备 
择 模 型 ， 并 根据 这 些 信息 选择 具体 的 分 析 方 法 。 
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1 引言 检验 而 言 ， 零 假设 的 形式 通常 为 模型 拟 合 良好 或 
模型 间 无 差异 ， 而 备 择 假设 则 通常 为 模型 拟 合 不 
结构 方程 模型 (Structural Equation Model, 良 或 模型 间 存 在 差异 。 针 对 此 类 假设 ,本 类 错误 
SEM) 是 心理 学 、 管 理学 、 社 会 学 等 学 科 中 重要 的 意味 着 研究 者 可 能 将 拟 合 不 良 的 模型 误 认 为 拟 合 
统计 工具 之 一 ,可 以 用 于 验证 性 因子 分 析 、 共 同 良好 的 模型 进行 结果 报告 , 王 阳 等 人 (2020) 则 指出 ， 
方法 偏差 检验 、 中 介 / 调 节 效 应 分 析 、 交 叉 灌 后 分 若 模型 设 定 有 误 ， 建 立 在 此 模型 基础 上 的 参数 估 
析 等 多 种 场景 。 有 研究 者 指出 ， 相 较 于 一 般 的 回 计 结 果 都 是 不 可 靠 的 。 在 这 种 情况 下 ， 开 类 错误 
归 分 析 , 结构 方程 模型 能 够 更 好 地 控制 测量 误差 ， 可 能 带 来 比 工 类 错误 更 大 的 危害 。 
也 支持 构建 复杂 的 多 变量 模型 ( 王 阳 等 ,2020)。 近年 来 ,有 一 些 研究 者 开始 逐渐 意识 到 使 用 
然而 , 在 使 用 结构 方程 模型 进行 分 析 时 ， 却 较 少 结构 方程 模型 作为 统计 工具 时 统计 检验 力 分 析 的 
有 研究 者 报告 统计 检验 力 的 相关 内 容 。 事 实 上 ， ”重要 性 ,并 开始 在 研究 中 报告 其 先 验 或 后 验 检验 
统计 检验 力 的 大 小 在 零 假设 显著 性 检验 (Null Fi, 或 者 至 少 在 讨论 部 分 指出 该 研究 可 能 在 统计 检 
Hypothesis Significance Testing, NHST) 中 有 着 重 验 力 上 存 有 不 足 (Mullen & Crowe, 2017; Hollerbach 
要 的 参考 价值 。 et al., 2018; Kornadt et al., 2019; Zhang & Zaman, 
统计 检验 力 指 的 是 某 假设 检验 能 够 正确 地 和 2020; Zhai et al., 2021)。 但 值得 注意 的 是 ， 目 前 仍 
绝 错误 零 假 设 的 概率 。 对 统计 检验 力 进 行 分 析 ， 有 许多 研究 者 在 使 用 结构 方程 模型 时 并 未 提 及 统 
可 以 帮助 研究 者 确定 合理 样本 量 ， 避 免 造成 大 量 计 检 验 力 的 有 关内 容 。 笔者 曾 以 全 文 包 含 “结构 方 
的 人 力 和 物力 的 浪费 ( 吴 艳 , EIB, 2011; 赵 礼 ， 程 模型 或 “SEM 模型 "为 条 件 ， 搜索 了 国内 心理 
ERE, 2019)。 此 外 ,对 评估 结构 方程 模型 拟 合 的 9? 学 顶级 刊物 《心理 学 报 》 中 近 5 年 (2017~2021) 内 
发 表 的 论文 ， 共 检索 出 51 篇 论文 ， 排 除 其 中 6 篇 
方法 类 研究 和 2 篇 元 分 析 研 究 , 余下 43 篇 中 提 及 
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和 具体 分 析 流 程 ， 并 进行 了 实例 演示 ， 最 后 对 相 
关 问 题 进 行 了 讨论 。 


2 结构 方程 模型 中 的 假设 检验 


假设 检验 的 技术 路 线 直 接 决 定 了 检验 力 分 析 
所 能 选取 的 技术 路 线 ， 结 构 方程 模型 中 假设 检验 
的 技术 路 线 大 致 可 分 为 两 种 : 基于 x? 分布 的 方法 
以 及 基于 模拟 (simulation) 或 重 抽样 (resampling) 技 
术 的 方法 。 
2.1 基于 * 分 布 的 结构 方程 模型 检验 

基于 分布 的 结构 方程 模型 检验 主要 包括 用 
于 评价 模型 拟 合 (或 肉 套 模型 比较 ) 的 x 检验、Ax 
检验 、 等 效 性 检验 、 基 于 模型 拟 合 指数 的 检验 以 
及 用 于 评价 模型 参数 (或 参数 的 函数 ， 如 中 介 效 应 ) 
的 Wald 检验 、 拉 格 朗 日 乘 子 检验 (Lagrange Multiplier 
test, LM 检验 ) 和 似 然 比 检验 (Satorra, 1989; Chou & 
Bentler, 1990; Gonzalez & Griffin, 2001; Yuan & 
Chan, 2016; Yuan et al., 2016; 王 阳 等 , 2020). 这些 
检验 涵盖 了 结构 方程 模型 相关 假设 检验 的 方 方 面 
面 , 其 至 可 以 说 , 目前 与 结构 方程 模型 有 关 的 参 
数 检验 方法 要 么 本 身 就 是 x 检验 , 要 么 可 以 通过 
适当 变形 转化 为 x 检验。 根据 构建 统计 量 时 的 细 
微 差别 ， 上 述 检 验 又 可 被 划分 为 三 个 子 类 。 
首先 是 x 检验 、Ax 检验、 等 效 性 检验 和 似 
然 比 检验 , 这 四 种 检验 都 以 结构 方程 模型 的 拟 合 
函数 为 基础 。 结 构 方程 模型 的 拟 合 函 数 衡 量 了 样 
本 协 方差 矩阵 S 与 假设 模型 对 应 的 协 方差 矩阵 Eo 
之 间 的 差距 ， 亦 即 变量 间 关 系 在 假设 模型 与 实际 
数据 中 的 差距 ( 匡 明 隆 ，2010; 王 阳 等 ,2020)。 
Muthén(2004) 指 出 ， 当 样本 来 自 假 设 模 型 所 对 应 
的 正 态 总 体 时 ， 拟 合 函 数值 的 N-I AEN 为 样本 量 ) 
随 着 N 的 增 大 渐 近 服从 中 心 化 好 分 布 。 另 有 学 者 
则 指出 ， 当 假定 样本 来 自 与 假设 模型 差距 小 于 给 
定 值 ( 记 作 0) 的 正 态 总 体 时 ， 此 时 拟 合 函数 值 的 
N-1 倍 渐进 服从 非 中 心 参数 为 (VW-7)Fo 的 非 中 心 
x 分 布 (Satorra & Saris, 1985; MacCallum et al., 
1996; Yuan & Chan, 2016)。 据 此 ,可 以 构建 出 零 假 
设 成 立时 x 检验 和 等 效 性 检验 所 对 应 的 统计 量 分 
布 。 而 AX 检验 和 似 然 比 检验 则 是 在 前 者 的 基础 
E, 增加 了 对 假设 模型 的 约束 条 件 , 通过 比较 添 
加 约束 条 件 前 后 实际 数据 与 假设 模型 的 差距 是 否 
显著 增 大 来 判定 添加 的 约束 条 件 是 否 “ 合 理 ”。 

第 二 类 是 基于 模型 拟 合 指数 的 检验 。 从 原理 


EVE, 这 类 检验 与 第 一 类 检验 差别 不 大 ， 其 本 质 
是 通过 拟 合 指数 的 定义 式 将 拟 合 指数 换算 为 非 中 
心 参数 ， 因 此 该 方法 仅 能 用 于 根据 非 中 心 参数 定 
义 的 拟 合 指数 (如 RMSEA, 、CFI、MFI、GFI 等 )。 
此 类 检验 有 两 种 相似 但 略 有 不 同 的 技术 路 线 ， 一 
种 将 基于 模型 拟 合 指数 的 零 假 设 换算 为 基于 非 中 
心 参数 的 零 假 设 ,， 选择 非 中 心 参数 的 临界 值 作为 
构建 x? 统 计量 分 布 的 依据 , MacCallum 等 人 (1996) 
和 MacCallum 等 人 (2006) 就 采用 了 这 种 做 法 。 另 
一 种 则 根据 Venables (1975) 的 研究 ， 以 样本 计算 
出 的 非 中 心 参数 估计 值 为 依据 构建 总 体 非 中 心 
参数 的 1 - wx 置信 区 间 ， 再 根据 定义 式 将 其 换算 
回 拟 合 指数 ， 从 而 得 到 拟 合 指数 的 1 - a 置信 区 
间 ， 进 而 完成 统计 推断 。 基 于 CFI 与 RMSEA 的 
等 效 性 检验 则 采用 了 此 类 做 法 (Yuan et al., 2016; 
Marcoulides & Yuan, 2017; 王 阳 等 , 2020)。 
第 三 类 是 Wald 检验 和 LM 检验 。 与 前 两 类 方 
法 不 同 的 是 ， 此 类 方法 不 再 关注 样本 的 拟 合 函 数 
值 本 身 ( 非 中 心 参数 也 是 关于 拟 合 函数 的 函数 ), 
而 是 通过 对 关键 参数 及 其 标准 误 的 估计 , 构建 z 
统计 量 进行 检验 。 其 中 Wald 检验 以 备 择 假 设 成 立 
时 的 模型 为 基础 完成 对 目标 参数 的 估计 和 检验 。 
当 处 理 参 数 函 数 形式 的 假设 时 ，Wald 检验 对 参数 
函数 在 其 真 值 邻 域内 进行 一 阶 Taylor 展开 求 得 其 
线性 近似 值 ， 进 而 完成 对 待 检验 指标 的 标准 误 估 
计 (Chou & Bentler, 1990; Gonzalez & Griffin, 
2001)。 此 时 Wald 检验 亦 称 为 Delta 法 ， 叶 宝 娟 和 
温 忠 麟 (2012) 兽 介绍 过 如 何 利用 该 方法 对 多 维 测 
验 合 成 信和 度 的 置信 区 间 进 行 估计 。 与 Wald 检验 相 
对 , LM 检验 以 零 假 设 成 立时 的 模型 为 基础 ， 将 假 
设 条 件 重 写 为 约束 函数 的 形式 ,并 构建 拉 格 朗 日 
函数 将 原 假设 检验 问题 变 为 检验 拉 格 朗 日 乘 数 是 
TA 0 的 问题 (Chou & Bentler, 1990)。 若 拉 格 朗 日 
乘 数 为 0, 则 意味 着 新 假设 条 件 的 加 入 没有 改变 
原 假设 模型 。 
22 ”基于 模拟 或 重 抽样 技术 的 结构 方程 模型 检验 
近年 来 ， 随 着 计算 机 技术 的 发 展 ,模拟 技术 
和 重 抽样 技术 越 来 越 多 地 出 现在 统计 学 领域 , 这 


' 本 文中 出 现 的 “GFI” 和 “AGFI” 指 的 分 别 是 “Unbiased 
goodness of fit index” 和 “Unbiased adjusted goodness of fit 
index”, {fiidE“Goodness of fit index” #il“Adjusted goodness 
of fit index”， 详 见 网 络 版 补充 材料 。 
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些 技术 被 广泛 地 运用 在 包括 参数 估计 与 假设 检验 
在 内 的 多 种 场合 。 不 同 于 基于 理论 分 布 的 经 典 方 
法 ， 此 类 方法 通常 以 模拟 或 重 抽 样 的 方式 得 到 目 
标 统计 量 的 经 验 分 布 , 并 以 经 验 分 布 代替 理论 
Ay Ap EAT Set HET, ALU EEA F H oe BE et 
的 理论 分 布 不 明确 或 难于 计算 的 情形 。 此 类 方法 
中 最 为 重要 的 两 种 分 别 是 基于 模拟 技术 的 Monte 
Carlo (MC) 方 法 ?与 基于 重 抽样 技术 的 Bootstrap 
方法 。 

MC 方法 最 早 由 Stanislaw Ulam 等 人 于 20 世 
纪 40 年 代 中 期 提出 ， 其 核心 思想 为 采用 模拟 的 方 
式 “ 重 现 ” 某 一 概率 模型 所 描绘 的 过 程 ， 以 模拟 得 
到 的 结果 作为 目标 问题 的 近似 解 。MacKinnon 等 
人 (2004) 在 评估 中 介 效 应 检验 的 M 法 和 经 验 M 法 
时 同时 引入 了 MC 法 。Preacher 和 Selig (2012) 进 
一 步 指出 MC 法 不 仅 可 以 用 于 线性 回归 和 路 径 分 
析 模 型 ， 也 可 直接 扩展 至 潜 变量 模型 。 该 方法 假 
定 模型 参数 服从 联合 正 态 分布 ， 其 均值 向 量 与 协 
方差 矩阵 通常 由 对 原始 模型 的 极 大 似 然 估计 确 
Eo MC 法 将 模型 参数 视 为 从 确定 的 联合 正 态 分 
布 中 抽样 的 过 程 ， 以 此 为 基础 进行 随机 抽样 ， 最 
终 基 于 抽样 结果 构建 中 介 效 应 的 经 验 分 布 并 进行 
统计 推断 (Preacher & Selig, 2012)。 此 外 , 在 结构 
方程 建 模 中 还 有 一 类 Monte Carlo 方法 一 一 马尔 
可 夫 链 蒙特 卡 洛 (Markov chain Monte Carlo, 
MCMC) 方 法 ， 此 类 方法 主要 用 于 贝 叶 斯 结构 方 
程 建 模 中 。 由 于 对 “参数 "这 一 重要 概念 的 基础 假 
设 不 同 ， 贝 叶 斯 统计 中 假设 检验 和 统计 检验 力 分 
析 的 过 程 均 与 基于 频率 学 派 的 传统 统计 学 有 着 较 
大 的 区 别 , 本 文 的 内 容 仍然 在 频率 学 派 的 框架 内 ， 
故 对 该 方法 不 做 过 多 介绍 。 

注意 到 ,基于 模拟 技术 的 结构 方程 模型 检验 
仍然 对 分 布 有 着 较 强 的 假设 , MC 法 假定 路 径 系数 
服从 联合 正 态 分 布 ,基于 贝 叶 斯 统计 的 MCMC 法 
需要 给 出 模型 参数 的 先 验 分 布 与 给 定 参 数 时 样本 
的 条 件 分 布 。 相 比 之 下 ,基于 重 抽 样 技术 的 检验 
对 样本 分 布 的 假设 更 弱 ，Efron (1979) 提 出 的 


”Monte Carlo 方法 实际 上 代表 的 是 以 计算 机 模拟 为 核心 的 
一 组 方法 , 为 了 和 后 文 统计 检验 力 分 析 的 Monte Carlo 方 法 
加 以 区 分 , 假设 检验 的 Monte Carlo 方法 本 文 统一 简称 为 
MC 方法 /MC 法 ,特此 说 明 。 


Bootstrap 法 ? 仅 要 求 样 本 是 从 某 一 确定 总 体 中 随 
机 抽样 得 到 的 ， 对 其 分 布 形态 并 无 更 多 要 求 。 不 
仅 如 此 , Bootstrap 法 同样 基于 “利用 经 验 分 布 代替 
理论 分 布 进行 统计 推 新 ”的 思路 , 在 面 对 一 些 理 
论 分 布 不 明 或 难以 计算 的 场合 表现 良好 。 由 于 
Bootstrap 法 具有 前 述 技术 优势 ， 近 年 来 ,该 方法 
已 逐渐 成 为 结构 方程 模型 假设 检验 问题 中 的 显 学 ， 
被 广泛 应 用 于 多 种 不 同 场合 。Zhang 和 Savalei 
(2016) 指 出 可 以 用 Bootstrap 法 估计 结构 方程 模型 
拟 合 指数 的 置信 区 间 。 Rosseel 在 其 编写 的 R 程序 
包 lavaan 中 给 出 了 求 结构 方 程 模 型 路 径 系 数 的 
Bootstrap 标准 误 的 程序 (Rosseel, 2021)。 更 有 研究 
者 利用 Bootstrap 法 估计 一 些 根据 结构 方程 模型 路 
径 系数 计算 得 到 的 更 复杂 的 指标 。 层 金 路 等 人 
(2005) 利 用 Bootstrap 法 对 问卷 合成 信和 度 的 置信 区 
间 进 行 了 估计 ; 张 涵 和 康 飞 (2016) 介 绍 了 基于 
Bootstrap 的 多 重 中 介 效 应 分 析 。 方 态 等 人 (2011) 
的 研究 及 温 忠 腾 和 叶 宝 娟 (2014a, 2014b) 的 研究 则 
分 别 介绍 了 Bootstrap 法 在 中 介 效 应 与 有 调节 的 中 
介 效 应 检验 中 的 应 用 。 
表 1 中 按照 不 同类 别 总 结 了 结构 方程 模型 涉 
及 到 的 常见 检验 的 应 用 场景 、 零 假设 、 统 计量 构 
建 方式 、 零 假设 成 立时 统计 量 对 应 的 分 布 以 及 使 
用 中 的 优 缺点 。 


3 结构 方程 模型 假设 检验 统计 检验 力 
的 分 析 方 法 


统计 检验 力 的 分 析 方法 与 分 析 对 象 所 采用 的 
技术 路 径 息 息 相 关 ， 结 构 方程 模型 假设 检验 统计 
检验 力 的 分 析 方 法 也 可 大 致 分 为 两 条 不 同 的 技术 
路 线 : 基于 x 分布 的 方法 和 基于 模拟 的 方法 。 其 
中 ,基于 从 分 布 的 方法 可 以 用 于 解决 使 用 经 典 方 
法 进行 结构 方程 模型 假设 检验 时 的 统计 检验 力 分 
析 问 题 ， 而 基于 模拟 的 方法 可 以 用 于 解决 一 般 的 
结构 方程 模型 假设 检验 的 统计 检验 力 分 析 问 题 。 
3.1 基于 分布 的 统计 检验 力 分 析 

顾名思义 ,基于 x 分布 的 方法 以 x0 分布 的 相 
关 理 论 为 基础 。 此 类 方法 可 以 追溯 到 Satorra 和 


3 事实 上 ，Bootstrap 法 还 分 为 参数 Bootstrap 法 和 非 参数 
Bootstrap 法 ,参数 Bootstrap 法 由 于 需要 了 解 随 机 变量 的 分 
布 型 态 , 在 实际 中 使 用 不 多 。 本 文中 关 不 加 说 明 , Bootstrap 
法 特 指 非 参数 Bootstrap 法 。 
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表 1 结构 方程 模型 中 的 假设 检验 
零 假设 成 立 
今 验 方 》 下 的 类 H 9 1 出 ak 
检验 方法 所 处 的 类 别 rela ne Hy 统计 量 构成 。 时 统计 量 对 。 。 优点 。 缺点 
oa me 应 的 分 布 
基于 x 入 检验 x 检 模型 拟 合 的 x 检验 模 Furo=0 (n-1):FuL 中 心 化 2 4) 计算 简便 ; 对 样本 
分 布 的 验 ”检验 型 拟 合 25E 布 工具 成 熟 。 ”的 分 布 
方 法 基于 刀 统 计量 检验 模 Fano >e (DFu Pè 好 分 有 较 强 
Sate 的 等 效 性 检验 型 拟 合 与 zo 间 的 差距 ti 的 假 
2 ha s 。 z = = z i 2 统计 量 
AX HRS EEE BEE Furos 一 Fyro4=0 (n1) (Fus 一 中心 化 x 分 的 构成 
检验 ”的 A 检验 mpg Toa = E Fura) 布 ee 
— KE 
模型 参数 的 似 检验 特 9 (0) = 0 (n1) (Fu se 中 心 化 X SP 限 制 ; 
然 比 检验 定 参 数 Lown =Xo nse  — Fur sax) 布 面 对 较 
ME of Z 检 Wald 检 验 ”检验 特 g(9)=0 sg(9) 标准 正 态 分 为 复杂 
分 布 理 验 定 参 数 Donar = Xo 无 约束 SE ,0) 布 的 参数 
论 的 检 KO M 
验 LM 检验 检验 特 和 =0 ES 标准 正 态 分 (如 中 
定 参 数 Dosar = Xo zak SE, 布 介 效 
冯 分 拟 合 指数 的 检 检验 模 5=0 给 定 的 拟 合 不 明 , 但 临界 
布 非 验 型 拟 合 了 = 指数 值 与 中 心 化 好 
中 心 分 布 有 关 
参 数 基于 拟 合 指数 检验 模 >d 给 定 的 拟 合 不 明 , 但 临界 
检验 ”的 等 效 性 检验 型 拟 合 ES 世间 的 差距 指数 值 与 非 中 心 闪 
大 于 /5x) 分 布 有 关 
基于 模 基于 模 Monte MC 法 检验 模 Fuo = 0 给 定 的 拟 合 经 验 分 布 。 “放宽 了 对 样 计算 较 
拟 或 重 拟 技 术 ae 型 拟 合 E= 指数 (包括 多 本 的 假设 ; 为 困 
由 样 技 的 方法 统计 量 ) 不 需 讨 论 目 难 、 计 
T 检验 特 s(9)=0 (0) Bei OR a Se 
法 定 参 数 Lown = Eo 无 约束 理论 分 布 
基于 重 Boots Bootstrap 检验 模 Fano =0 给 定 的 拟 合 经 验 分 布 灵活 度 高 ， 法 要求 
抽样 技 型 拟 合 = Zo 指数 (包括 k 对 统计 量 的 使 用 者 
术 的 方 统计 量 ) 构造 几乎 没 有 编程 
i 检验 特 (0) = 0 (0) 经 验 分 布 。 ”有 限制 ; 对 技术 
定 参 数 Lowe = Eo xa 样本 分 布 要 
求 宽松 


Saris 早年 的 一 系列 研究 (Satorra & Saris, 1985; 
Satorra, 1989; Satorra et al., 1991)。Satorra (1989) 
指出 ， 当 前 述 假设 检验 的 加 不 成 立时 ,其 对 应 统 


计量 (对 Wald 检验 和 LM 检验 而 言 是 统计 量 的 平 
方 ) 服 从 非 中 心 阅 分布, 且 其 相应 的 非 中 心 参数 可 


以 通过 总 体 协 方差 矩阵 真 值 拟 合 零 假 设 模型 计算 


出 的 统计 量 得 到 。 进 而 可 以 通过 计算 备 择 假设 成 
SEARLE PEA Serbs At 


E 绝 域 的 概率 求 得 给 定 


检验 的 统计 检验 力 。 设 零 假设 和 备 择 假设 成 立时 
对 应 的 非 中 心 参数 分 别 为 8。 和 61， 则 统计 检验 力 
1 一 B 由 下 式 给 出 : 


1- p= Pirin < Tuc} 


Bu - B= PL 73,5 > 
(对 应 除 等 效 性 检验 外 的 
HREA df, 非 中 心 参数 为 5 的 非 


其 中 72, RR A 
中 心 卡 方 分 布 , Pwzc 由 下 式 确定 : 


Ps > Tuc} =g 


(拒绝 域 在 左 


I, 对 应 等 效 性 检验 ) 


(拒绝 域 在 左 侧 , 对 应 等 效 性 检验 ) (La) 
Tyc 拒绝 域 在 右 侧 ， 
其 他 检验 ) (1b) 


(2a) 


BLP Zio, S Truc} = eGR REM, 


对 应 除 等 效 性 检验 外 的 其 他 检验 ) 


I 


(2b) 
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由 于 56 在 可 确定 时 已 经 可 以 算出 , 因此， 基 
于 分 布 的 统计 检验 力 分 析 的 关键 在 于 估计 61。 
根据 估计 非 中 心 参数 的 方法 不 同 ， 此 类 方法 又 可 
分 为 利用 总 体 协 方差 矩阵 E 进行 估计 的 Satorra- 
Saris 法 和 利用 拟 合 指数 (如 RMSEA CFI, MFI 
等 ) 进 行 估 计 的 MacCallum 法 。 
3.1.1 Satorra-Saris 法 

Satorra-Saris 法 主要 基于 Satorra 和 Saris 先前 
的 研究 ， 其 具体 步骤 如 下 (Satorra & Saris, 1985; 
Satorra et al., 1991): 

第 一 步 : 定义 备 择 假设 成 立 情况 下 的 模型 ， 
指定 备 择 假设 成 立 情况 下 的 理论 总 体 协 方差 矩 
4X, 

第 二 步 : 设 定 样本 容量 为 真实 样本 容量 N, 
利用 三 拟 合 零 假设 模型 。 

第 三 步 : 将 第 二 步 中 计算 得 到 的 统计 量 (针对 
六 分 布 的 统计 量 为 统计 量 本 身 , 针对 标准 正 态 分 
布 的 统计 量 为 统计 量 的 平方 ) 作 为 61 的 估计 ,代入 
(1a) 或 (1b) 完 成 统计 检验 力 计算 。 

使 用 Satorra-Saris 法 估计 达到 某 统计 检验 力 
所 需 样本 量 的 过 程 与 估计 统计 检验 力 类 似 ， 只 
是 第 二 步 中 需要 反复 设 定 样本 容量 N 进行 多 次 
估计 ， 并 计算 不 同样 本 容量 对 应 的 统计 检验 力 ， 
最 终 从 中 选取 满足 目标 统计 检验 力 所 需 的 最 小 
样本 容量 。 

Satorra-Saris 法 通过 定义 不 同 于 零 假 设 模型 
的 备 择 模 型 重 构 总 体 协 方差 矩阵 ， 并 据 此 计算 统 
计 检 验 力 。 其 中 ,对 备 择 模型 进行 定义 的 过 程 不 
仅 需要 研究 者 给 出 明确 的 模型 结构 ， 还 需要 研究 
者 指定 大 量 的 参数 (如 路 径 系 数 、 测 量 残 差 等 )。 因 
此 ， 从 应 用 的 角度 考虑 ，Satorra-Saris 法 更 适合 在 
研究 者 对 备 择 模型 的 了 解 程度 较 深 的 情况 下 使 用 
(如 针对 模型 参数 的 检验 )。 

3.1.2 MacCallum 法 

Satorra-Saris 法 对 备 择 模型 过 强 的 假设 从 客 
观 上 提高 了 该 方法 的 使 用 难度 ,也 限制 了 其 适用 
范围 。MacCallum 等 人 (1996) 给 出 了 另外 一 种 估计 
韭 中 心 参数 的 方法 , 该 方法 同样 基于 Satorra 和 
Saris (1985) 给 出 的 公式 ， 但 在 估计 Furo 时 ， 
MacCallum 等 人 给 出 了 一 个 基于 RMSEA 的 估计 
方法 。 他 们 指出 ,使 用 样本 拟 合 函 数 估计 值 Fu 
作为 对 总 体 参数 Po 的 估计 是 有 偏 的 ， 且 Furo- 
Fur Z BA Wh F (MacCallum et al., 1996): 


Furo = Fy = of df -RMSEA’ (3) 


-1 

据 此 , MacCallum 等 人 提出 基于 RMSEA 指标 
的 统计 检验 力 分 析 步 又 (MacCallum et al., 1996): 

第 一 步 : 根据 零 假 设 和 备 择 假设 成 立时 对 应 
的 RMSEA 值 、 模 型 自由 度 以 及 样本 量 分 别 计算 
出 零 假设 和 备 择 假设 成 立时 的 非 中 心 参数 Sy 和 61。 

第 二 步 : 代入 (1a) 或 (1b) 完 成 统计 检验 力 计算 。 

与 Satorra-Saris 法 类 似 ， 当 使 用 MacCallum 
法 估计 达到 某 统计 检验 力 所 需 样本 量 时 ， 也 只 需 
要 在 第 一 步 中 反复 设 定 不 同 的 N 对 统计 检验 力 进 
行 多 次 估计 ,最终 从 中 选取 满足 目标 统计 检验 力 
所 需 的 最 小 样本 容量 。 

Kim (2005) 对 MacCallum 法 进行 了 进一步 简 
化 和 扩展 ,给 出 了 一 个 快速 求 取 给 定 统计 检验 力 
情况 下 所 需 样本 量 的 公式 : 
CE Oe ey 
df -RMSEA? 

其 中 qf 为 模型 自由 度 , RMSEA 为 备 择 假设 成 
立时 对 应 的 RMSEA 值 , 5 为 显著 性 水 平 为 a、 统 
计 检 验 力 为 1-B 情况 下 对 应 的 非 中 心 参数 。 非 中 
心 参数 可 以 通过 查 表 或 者 使 用 统计 软件 计算 得 到 ， 
Kim (2005) 在 他 们 的 文章 最 后 同时 给 出 了 Ho: 
Pixzo=0 时 利用 计算 机 快速 求 取 8 的 SPSS 与 SAS 
脚本 。 此 外 , Kim (2005) 在 他 的 文章 中 还 提出 了 基 
F CFI, MFI 与 GFI 的 计算 公式 。 

相 比 于 Satorra-Saris 法 , MacCallum 法 采用 拟 
合 指 数 对 模型 的 拟 合 情况 进行 了 某 种 意义 上 的 
“打包 ”， 即 利用 拟 合 指数 对 假设 模型 与 实际 数据 
的 差距 进行 了 概括 ， 使 研究 者 免 于 设 定 复杂 的 备 
择 模型 ， 在 分 析 模 型 拟 合 检验 的 统计 检验 力 时 
MacCallum 法 仅 需要 模型 自由 度 、 零 假设 和 备 择 
假设 成 立时 对 应 的 拟 合 指数 值 [ 若 对 关于 CFI 的 检 
验 进行 检验 力 分 析 时 还 需要 提供 基础 模型 (base 
model) MIU PBUH Fre 与 自由 度 dfg, 关于 GFI 
的 检验 还 需要 提供 协 方差 矩阵 的 秩 p]， 大 大 简化 
了 统计 检验 力 的 分 析 难 度 。 但 由 于 模型 参数 与 拟 
合 指数 之 间 存 在 复杂 的 非 线 性 关系 ，MacCallum 
法 无 法 分 析 模 型 参数 检验 的 统计 检验 力 。 
3.2 ”基于 模拟 的 统计 检验 力 分 析 

采用 模拟 方法 对 结构 方程 模型 假设 检验 进行 
统计 检验 力 分 析 早 年 可 追溯 到 Muthén 与 Muthén 
(2002) 的 研究 。 采 用 模拟 方法 进行 统计 检验 力 分 


(4) 


N RMSE. 
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析 的 思路 与 采用 模拟 方法 进行 结构 方程 模型 假 
设 检验 的 思路 相似 ， 即 利用 模拟 的 方式 确定 
P( 拒 绝 H610=0') 的 经 验 值 ,具体 而 言 ,此 类 方法 
通常 采用 模拟 的 方式 从 参数 为 0 的 特定 分 布 中 生 
成 多 组 随机 样本 ， 而 后 针对 每 组 随机 样本 执行 同 
样 的 假设 检验 程序 ， 最 后 统计 全 部 已 进行 的 假设 
检验 中 拒绝 零 假 设 的 频率 ， 以 此 作为 对 该 检验 统 
计 检 验 力 的 分 析 结 果 。Muthén 与 Muthén (2002) 
的 研究 给 出 了 利用 Monte Carlo 模拟 计算 基于 经 
典 方法 的 结构 方程 模型 假设 检验 的 统计 检验 力 的 
分 析 过 程 与 Mplus 程序 。Zhang (2014) 的 研究 指出 
Monte Carlo 模拟 法 同样 可 以 估计 基于 Bootstrap 
方法 的 中 介 效 应 检验 的 统计 检验 力 , 并 给 出 了 一 个 
基于 R 语言 bmem 程序 包 的 示例 程序 。Schoemann 
等 人 (2017) 的 研究 在 实现 对 中 介 模 型 的 统计 检验 
力 分 析 时 使 用 了 Monte Carlo 模拟 法 估计 基于 MC 
法 的 中 介 效 应 检验 的 统计 检验 力 。 总 的 来 说 , 使 
用 Monte Carlo 模拟 法 进行 统计 检验 力 分 析 的 步 
又 大 体 如 下 (Muthen & Muthén, 2002; Zhang, 2014; 
Schoemann et al., 2017): 

第 一 步 : 根据 备 择 假设 成 立时 的 情况 设 定 月 
于 生成 模拟 数据 的 均值 向 量 (模型 不 考虑 均值 结 
构 时 为 零 向 量 ) 与 协 方差 矩阵 。 此 步骤 中 既 可 以 
如 Schoemann 等 人 (2017) 直 接 指 定 协 方差 矩阵 ， 
也 可 如 Zhang (2014) 通 过 定义 模型 来 间接 指定 协 
方差 和 矩阵。 此外， 如 果 研 究 者 希望 模拟 非 正 态 
分 布 的 数据 ,通常 还 需要 给 出 目标 分 布 的 峰 度 
和 偏 度 。 

第 二 步 : 确定 需要 进行 Monte Carlo FAW AYE 
复 轮 数 R 与 样本 容量 N, Zhang (2014) 建 议 重复 轮 
数 尺 通常 不 应 少 于 1000, 

第 三 步 : 在 每 一 轮 Monte Carlo 模拟 中 ,重复 
执行 下 述 步 又: (1) 根 据 给 定 的 分 布 参数 与 样本 量 
生成 服从 特定 分 布 的 随机 数 作 为 模拟 样本 ; (2) 对 
该 模拟 样本 进行 结构 方程 建 模 并 进行 检验 ， 具体 


ag 


ran 


值 、 参 数 估 计 结 果 的 平均 值 和 标准 差 、 参 数 估计 
标准 误 的 平均 值 以 及 参数 95% 置 信 区 间 对 理论 真 
值 的 覆盖 率 ,Muthén 和 Muthén (2002) 进 一 步 提出 
衡量 Monte Carlo 模拟 结果 是 否 可 靠 的 标准 : (1) 模 
型 中 任何 参数 及 其 标准 误 真 值 与 估计 值 之 间 的 相 
对 误差 不 超过 10% (参数 估计 结果 的 标准 差 视 为 
标准 误 的 真 值 ); (2) 模 型 中 关键 参数 及 其 标准 误 真 
值 与 估计 值 之 间 的 相对 误差 不 超过 5%; (3) 参 数 的 
95% 置 信 区 间 对 理论 真 值 的 覆盖 率 应 当 在 0.91 至 
0.98 之 间 。 

当 使 用 Monte Carlo 法 估计 达到 某 统计 检验 
力 所 需 样本 量 时 , 也 只 需要 在 第 一 步 中 反复 设 定 
不 同 的 V 对 统计 检验 力 进 行 多 次 估计 ， 而 后 通过 
线性 插值 法 对 所 需 样本 量 进行 估计 。 与 基于 ?分 
布 的 统计 检验 力 分 析 方 法 相 比 ,结构 方程 模型 统 
计 检 验 力 的 Monte Carlo 分 析 由 于 不 需要 计算 目 
标 统计 量 的 理论 分 布 ， 因 此 可 以 广泛 适用 于 对 各 
种 检验 进行 统计 检验 力 分 析 的 情形 , 在 给 出 峰 度 
和 偏 度 时 ， 该 方法 还 可 以 在 一 定 程 度 上 分 析 总 体 
分 布 为 非 正 态 分 布 的 情形 。 但 结构 方程 模型 统计 
检验 力 分 析 的 Monte Carlo 法 需要 研究 者 至 少 给 
出 用 以 生成 模拟 数据 的 协 方差 矩阵 ， 这 意味 着 研 
究 者 需要 对 备 择 假 设 成 立 下 数据 的 分 布 情况 有 足 
够 的 认识 。 另 外 ,该 方法 本 身 对 研究 者 的 编程 水 
平 要 求 较 高 ， 且 计算 量 较 大 ， 计 算 起 来 较为 耗 时 ， 
在 一 定 程 度 上 影响 了 该 方法 的 使 用 和 推广 。 
3.3 ”结构 方程 模型 统计 检验 力 分 析 的 程序 资源 

为 了 方便 研究 者 对 统计 检验 力 进行 分 析 ， 多 
名 研究 者 提供 了 在 线 或 离线 的 相关 程序 ,基于 
分 布 的 统计 检验 力 分 析 方 法 中 , Preacher 和 Coffman 
(2006) 提 供 了 一 个 基于 MacCallum 法 的 统计 检验 
力 在 线 计算 程序 ， 可 以 计算 给 定 样本 量 情况 下 ， 
单个 模型 拟 合 或 藤 套 模型 拟 合 差异 的 统计 效力 或 
达到 某 统计 效力 所 需 的 样本 量 。Moshagen 提供 了 


的 检验 方法 可 以 是 表 1 中 的 任意 一 种 。 

第 四 步 : 统计 出 每 一 轮 统 计 检 验 的 结果 ， 设 
在 R 轮 检 验 中 有 7 轮 结果 拒绝 了 原 假设 ， 则 最 终 
的 统计 检验 力 为 二 。 


此 外 ， 为 了 说 明 Monte Carlo 模拟 结果 的 准确 
性 , Muthén 和 Muthén (2002) 建 议 研究 者 在 报告 结 
果 的 同时 报告 结构 方程 模型 中 各 参数 的 理论 真 


个 R 程序 包 semPower, 该 程序 包 基 于 Moshagen 
和 Erdfelder (2016) 的 研究 ， 可 以 计算 基于 Furos 

RMSEA、McDonald 拟 合 指数 、GFI 与 AGFI 的 统 
计 检 验 力 ， 当 指定 零 假 设 和 备 择 假设 成 立 情况 下 
对 应 的 协 方差 矩阵 时 ， 该 程序 包 也 可 使 用 Satorra- 
Saris 法 分 析 统 计 检 验 力 。Jak 等 人 (2021) 给 出 
了 一 个 基于 Shiny 开发 的 在 线程 序 power4SEM 
(网 址 为 https://sjak.shinyapps.io/power4SEM/), 
该 程序 同时 支持 Satorra-Saris 法 和 MacCallum 法 ， 


A 
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Jak 等 人 在 他 们 的 论文 中 详细 介绍 了 该 程序 的 使 用 
方式 。 

基于 模拟 的 统计 检验 力 分 析 方 法 中 , Muthén 
All Muthén (2002) 提 供 了 使 用 Monte Carlo 法 分 析 
结构 方程 模型 统计 检验 力 的 Mplus 程序 框架 。 
Zhang 和 Wang (2020) 提 供 了 可 以 用 来 分 析 中 介 效 
应 统计 检验 力 的 R 语言 程序 包 bmem。 该 程序 包 
使 用 Monte Carlo 法 进行 统计 检验 力 分析 ， 可 以 对 
中 介 效 应 的 Sobel 检验 与 Bootstrap 检验 的 统计 检 
验 力 进行 分 析 。 特 别 的 是 ， 该 程序 包 还 提供 给 研 
究 者 自主 设 定 分 布 的 峰 度 与 偏 度 的 功能 以 模拟 总 
体 为 非 正 态 分 布 时 的 情形 。Schoemann 等 人 (2017) 
给 出 了 一 个 使 用 Monte Carlo 法 分 析 中 介 效 应 统 
计 检 验 力 的 在 线程 序 ( 网 址 为 https://schoemanna. 


的 平行 中 介 模 型 和 两 个 中 介 变 量 的 链 式 中 介 模 型 
三 种 情况 下 的 中 介 效 应 统计 检验 力 分 析 。Wang 和 
Rhemtulla (2021) 则 给 出 了 一 个 分 析 结 构 方程 模型 
参数 检验 统计 检验 力 的 在 线程 序 (网 址 为 : https:// 
yilinandrewang.shinyapps.io/pwrSEM/)_ 

表 2 中 总 结 了 结构 方程 模型 统计 检验 力 分 析 
方法 的 适用 场合 、 需 要 的 计算 量 以 及 相应 的 程序 
资源 。 


4 应 用 实例 


下 面 用 一 个 例子 演示 在 实际 研究 中 如 何 灵活 
使 用 各 种 方法 对 结构 方程 模型 的 统计 检验 力 进 行 
分 析 。 注 意 到 , 统计 检验 力 的 分 析 包 括 先 验 分 析 
和 后 验 分 析 ， 二 者 对 应 研究 过 程 中 的 不 同时 期 ， 
且 不 能 相互 蔡 代 ( 吴 艳 , TLE BR, 2011; 赵 礼 , E 
晖 ， 2019)。 但 先 验 分 析 与 后 验 分 析 在 分 析 步 又 上 


差异 不 大 ， 且 亦 有 人 研究 者 指出 , 在 NHST 中 , 统计 
检验 力 是 p 值 的 函数 , 在 已 知 p 值 的 情况 下 统计 
检验 力 没有 给 出 额外 信息 ， 因 此 建议 无 需 进 行 后 
验 分 析 (Hoenig & Heisey, 2001)。 考虑 到 篇 幅 所 限 ， 
本 文 仅 展示 先 验 分 析 的 过 程 ， 后 验 分 析 的 部 分 参 
见 网 络 版 补充 材料 。 演 示 的 样 例 来 自 Zhai 等 人 
(2021) 的 研究 ,该 研究 调查 了 新 冠 疫情 期 间 个 体 
的 心理 健康 与 创伤 后 成 长 的 基本 情况 ， 并 讨论 了 
在 控制 疫情 相关 压力 和 年 龄 的 情况 下 ,情绪 创造 
力 对 个 体 心理 健康 的 保护 作用 以 及 对 个 体 创 伤 后 
成 长 的 促进 作用 和 可 能 的 作用 机 制 。 具 体 的 模型 
设置 参见 图 1, 统计 软件 采用 R 3.6.3 的 lavaan 程 
序 包 和 semPower 程 序 包 ,本 部 分 全 部 及 代码 参见 
网 络 版 补充 材料 部 分 。 


T 


= shinyapps.io/me_power_med/), ELF AEF Shiny 从 检验 目的 的 角度 出 发 , 使 用 结构 方程 模型 
a FR, 可 用 来 处 理 简单 中 介 模 型 、 两 个 中 介 变 量 。 ”进行 验证 性 因子 分 析 的 主要 目的 是 通过 对 模型 拟 


合 的 评价 考察 样本 数据 是 否 符合 量 表 的 理论 结构 ， 
即 考察 量 表 题 项 与 测量 因子 间 的 关系 是 否 同 研究 
者 的 假设 所 一 致 ， 而 使 用 结构 方程 模型 进行 中 介 
效应 分 析 ， 特 别 是 对 显 变量 链 式 中 介 模 型 进行 中 
介 效 应 分 析 时 ,人 研究 者 的 目的 则 更 多 指向 对 路 径 
系数 函数 ( 即 中 介 效 应 ) 的 分 析 与 检验 。 此外， 注意 
到 图 1 中 的 模型 为 饱和 模型 有 研究 者 指出 ， 此 
时 模型 为 完美 拟 合 ， 可 不 关注 模型 拟 合 指数 ， 仅 
讨论 路 径 系 数 即 可 ( 张 莉 等 , 2019)。 综 合 考 虑 检 
验 目 的 与 实际 情况 ， 对 该 研究 的 统计 检验 力 进行 
分 析 时 ， 需 要 对 验证 性 因子 分 析 部 分 模型 拟 合 检 
验 的 统计 检验 力 与 链 式 中 介 模 型 部 分 模型 参数 检 
验 的 统计 检验 力 进 行 分 析 ， 链 式 中 介 效 应 部 分 选 
择 “ 情 绪 创 造 力 一 领悟 社会 支持 一 情绪 调节 自我 
效能 感 一 压力 后 成 长 "路径 的 中 介 效应 值 为 分 析 
对 象 。 


= 


表 2 结构 方程 模型 统计 检验 力 分 析 方法 


作 期 乔 


是 否 需要 给 出 对 应 假设 检验 
PEER 是 否 需 要 给 HIRSI pay 器 应 资 ， 
分 析 方 法 协 方差 矩阵 的 技术 路 线 假设 检验 的 目的 计算 量 程序 资源 
Satorra-Saris 法 是 仅 经 典 方法 均 可 , 但 更 适合 模型 参数 BY) power4SEM 
检验 或 嵌 套 模型 比较 we 
MacCallum 法 否 仅 经 典 方法 模型 拟 合 检验 最 小 power4SEM 
semPower 
Preacher 在 线程 序 
Monte Carlo 法 是 均 可 均 可 大 Bmem 


mc_power_med 
pwrSEM 


2124 心理 科学 进展 


ChinaXiv 合 作 期 干 


第 30 郑 


领悟 社会 支持 7 情绪 调节 自我 效能 感 


i 


图 1 样 例 的 结构 方程 模型 示意 图 
(虚线 表示 该 变量 为 控制 变量 ) 


41 ”模型 拟 合 检验 的 先 验 检验 力 分 析 
4.1.1 ”使 用 Satorra-Saris 法 分 析 模 型 拟 合 检验 的 
先 验 检验 力 

在 进行 先 验 检验 力 分 析 时 ， 首 先 要 设 定 研 究 
对 象 理论 上 的 总 体 协 方差 矩阵 (下 简称 理论 协 方 
差 和 矩阵)， 即 备 择 假设 成 立时 对 应 的 协 方差 矩阵 。 
与 上 检验 ,方差 分 析 的 先 验 检验 力 分 析 类 似 , 在 结 
构 方 程 模型 相关 检验 的 先 验 检验 力 分 析 中 ,理论 
协 方差 矩阵 也 可 以 通过 参考 前 期 理论 、 已 发 表 的 
类 似 成 果 或 研究 者 提出 的 标准 所 确定 ( 吴 艳 ， 温 忠 
W, 2011)。 通 常情 况 下 ,建议 研究 者 优先 参考 已 有 
类 似 研 究 中 报告 的 结果 设 定理 论 协 方差 矩阵 ， 当 
缺乏 既往 研究 或 既往 研究 信息 不 充足 时 , 研究 者 
可 以 参考 前 期 理论 或 其 他 研究 者 提出 的 标准 所 确 
定 。 此 外 , Jak 等 人 (2021) 建 议 研 究 者 在 设 定理 论 
协 方差 矩阵 时 ， 可 以 以 标准 化 的 形式 指定 参数 值 ， 
且 参 考 相关 研究 领域 中 * 小 、 中 、 大 ”所 对 应 的 
经 验 值 代入 待定 参数 。 同 时 ， 由 于 通常 情况 下 确 
定 协 方差 矩阵 所 需要 的 参数 数量 较为 庞大 , Jak 等 
人 (2021) 还 建议 研究 者 在 设 定理 论 协 方 差 矩 阵 时 
需要 考虑 那些 在 理论 上 能 够 影响 拒绝 或 接受 Ho 
的 重要 参数 ， 例 如 路 径 系 数 、 误 差 相 关 、 洪 变量 
相关 等 。 

通过 查询 资料 发 现 : 既往 研究 中 , 文书 锋 等 
人 (2009) 的 研究 较为 全 面 地 报告 了 可 以 用 于 生成 
理论 协 方差 矩阵 的 全 部 参数 ， 且 该 研究 样本 量 较 
大 (验证 性 因子 分 析 部 分 N = 818)。 因 此 , 在 之 后 
的 分 析 中 ， 以 该 模型 为 基准 生成 理论 协 方差 矩 
阵 。 本 例 中 控制 1 类 错误 水 平 a 为 0.05, 求 得 统 
计 检 验 力 为 0.95 时 对 应 的 样本 量 约 为 134 Ao HE 
外 ,对 不 熟悉 R 的 研究 者 ， 也 可 以 使 用 常用 的 结 
构 方 程 建 模 软件 (如 LISREL EQS, AMOS, Mplus 
等 ) 结 合 G*Power 进行 分 析 ， 有 具体 的 步骤 参见 网 络 
版 补充 材料 部 分 。 


ua 


4.1.2 ”使 用 MacCallum 法 分 析 模 型 拟 合 检验 的 
先 验 检验 力 

相 比 Satorra-Saris 法 , 使 用 MacCallum 法 分 
析 模 型 拟 合 检验 先 验 检验 力 的 过 程 要 简单 许多 ， 
只 需要 确定 零 假 设 和 备 择 假设 成 立 情况 下 对 应 的 
拟 合 指数 、 模 型 自由 度 df、 显 著 性 水 平 a 和 统计 
检验 力 1 一 B 五 个 指标 即 可 完成 分 析 , 在 设 定 备 择 
假设 对 应 的 拟 合 指数 时 ,研究 者 可 参考 既往 研究 
中 提出 的 关于 该 指数 “大 小 ”的 标准 。 例 如 对 
RMSEA 而 言 ， 有 研究 者 指出 ， RMSEA 小 于 0.05 
为 良好 拟 合 , 0.05 至 0.08 为 合理 拟 合 , 0.08 至 0.1 
为 可 接受 的 拟 合 , 大 于 0.1 为 拟 合 较 差 (MacCallum 
et al., 1996; Yuan & Chan, 2016; Jak et al., 2021)。 
对 于 本 文中 给 出 的 例子 ,控制 [类 错误 水 平 a 为 
0.05, AHA a BE RMSEA 的 真 值 为 0.08 
时 ,正确 拒绝 该 模型 的 几率 不 小 于 95%， 此 时 使 
用 semPower 包 提 供 的 semPower.aPriori 函数 可 以 
计算 得 出 需要 的 样本 量 为 134 人 。 
4.1.3 ”使 用 Monte Carlo 法 分 析 模 型 拟 合 检验 的 

先 验 检验 力 

使 用 Monte Carlo 法 分 析 模 型 拟 合 检验 先 验 
检验 力 时 确定 理论 协 方差 矩阵 的 步骤 和 原则 同 
Satorra-Saris 法 一 臻 ,这 里 不 再 歼 述 。 本 例 中 ,以 
等 效 性 检验 为 分 析 对 象 ， 控 制 工 类 错误 水 平 a 为 
0.05, 重复 次 数 R 为 1000。 分 别 取 N X 101, 201, 
301 的 情形 ,计算 在 给 定 N 的 情况 下 1000 次 重复 
试验 中 模型 拟 合 指数 落 在 以 RMSEA = 0.08 为 经 
验 临界 值 的 等 效 性 检验 拒绝 域 的 比率 。 需 要 注意 
的 是 ， 在 Monte Carlo 模拟 的 过 程 中 ， 有 些 样本 会 
导致 结构 方程 模型 在 计算 过 程 中 出 现 迭 代 不 收 
敛 、 残 差 方 差 为 负数 、 湾 变量 协 方差 矩阵 非 正 定 
等 “不 合理 ”的 情况 ， 即 出 现 了 不 恰当 解 ， 这些 解 
应 在 在 后 续 的 计算 过 程 中 予以 排除 。 本 例 中 , N 分 
别 为 101 和 201 时 对 应 的 恰当 解 比 例 分 别 为 
99.6% 和 100%， 对 应 的 统计 检验 力 分 别 为 0.724 
和 0.996; 对 应 的 各 参数 估计 情况 、 相 对 误差 以 及 
95% 置 信 区 间 履 盖 率 参见 网 络 版 补充 材料 。 本 例 
中 全 部 参数 均值 和 标准 误 的 估计 误差 均 不 超过 
10%, 95% 置 信 区 间 和 覆盖 率 均 处 于 0.91 至 0.98 之 间 ， 
提示 Monte Carlo 模拟 结果 的 可 信和 度 较 高 。 进 一 步 
根据 线性 插值 法 可 求 得 统计 检验 力 为 0.95 时 对 应 
的 样本 量 约 为 185 人 。 针 对 同一 个 假设 检验 进行 
统计 检验 力 分 析 时 , Monte Carlo 方法 得 到 的 结果 
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看 似 与 Satorra-Saris 法 有 较 大 差异 , 但 这 主 
由 于 线性 插值 带 来 的 误差 。 事实 上 ,车 通过 缩减 
采样 点 间距 的 方式 提高 线性 插值 的 精度 ， 可 重新 
求 得 待 估 样 本 量 约 为 151 Ao 
42 ”模型 参数 检验 的 先 验 检验 力 分 析 
4.2.1 ”使 用 Satorra-Saris 法 分 析 模 型 参数 检验 的 
先 验 检验 力 
使 用 Satorra-Saris 法 分 析 模 型 参数 检验 先 验 
检验 力 的 第 一 步 仍 然 是 设 定理 论 协 方差 矩阵 。 诚 
然 经 典 的 Satorra-Saris 法 强调 通过 备 择 模 型 生成 
理论 协 方差 矩阵 , 但 Satorra 和 Saris (1985) 给 出 的 
关于 该 方法 原理 的 讨论 中 并 没有 限制 理论 协 方差 矩 
阵 的 生成 模式 。 因 此 ,参照 Schoemann 等 人 (2017) 
的 方式 直接 设 定理 论 协 方差 矩阵 也 并 无 不 可 。 
由 于 既往 研究 中 缺乏 可 以 用 来 直接 重 构 理 论 
协 方差 矩阵 的 数据 ， 而 通过 元 分 析 技 术 虽 然 能 够 
在 一 定 程 度 上 基于 现 有 研究 构建 出 理论 协 方差 矩 
阵 , 但 这 种 方法 需要 耗费 较 多 的 人 力 物力 , 单纯 
用 来 进行 统计 检验 力 分 析 未 免 有 些 浪费 。 因 此 ， 
本 文 在 进行 理论 协 方差 矩阵 生成 时 ， 选 择 根据 相 
关 研 究 领 域 中 “小 、 中 、 大 ”所 对 应 的 经 验 值 进行 
设 定 。 此 外 ,由 于 尺度 变换 不 影响 对 本 例 中 所 关 
心 的 系数 ( 显 变量 间 的 路 径 系数 ) 显 著 性 的 判定 。 因 
此 , 在 指定 理论 协 方差 矩阵 时 可 进一步 将 其 简化 
为 指定 相关 矩阵。 参考 Cohen (1988) 的 观点 ， 相 关 
系数 在 0.1, 0.3, 0.5 时 分 别 对 应 “小 ”、“ 中 ”、“ 大 ” 
的 效应 量 。 在 本 例 中 , 假定 与 控制 变量 有 关 的 相 
关系 数 有 “小 ”的 效应 量 ， 关 注 的 变量 的 相关 系 
数 有 “中 ”的 效应 量 ， 以 此 设 定理 论 协 方差 矩阵 。 
以 对 “情绪 创造 力 一 领悟 社会 支持 一 情绪 调节 
自我 效能 感 一 压力 后 成 长 "这 一 路 径 中 介 效 应 的 
Wald 检验 为 分 析 对 象 ， 控制 工 类 错误 水 平 a 为 
0.05， 求 得 统计 检验 力 为 0.95 时 对 应 的 样本 量 约 
为 795 人 。 
4.2.2 ”使 用 Monte Carlo 法 分 析 模 型 参数 检验 的 

先 验 检验 力 

本 例 中 ,以 对 “情绪 创造 力 一 领悟 社会 支持 
一 情绪 调节 自我 效能 感 一 压力 后 成 长 "这 一 路 径 
中 介 效 应 的 wald 检 验 为 分 析 对 象 ， 控 制 1 类 错误 
水 平 a 为 0.05， 分 别 取 N X 301, 401, 501, 601 的 
情形 进行 模拟 , 得 到 的 结果 全 部 为 恰当 解 ， 对 应 
的 统计 检验 力 为 0.613，0.884，0.960，0.985， 对 应 
的 各 参数 估计 情况 、 相 对 误差 以 及 95% 置 信 区 间 


h 


覆盖 率 参 见 网 络 版 补充 材料 。 本 例 中 全 部 参数 的 
估计 误差 均 不 超过 10%， 关 键 参数 的 估计 误差 均 
不 超过 5%,95% 和 置信 区 间 和 覆盖 率 均 处 于 0.91 至 
0.98 之 间 ， 提示 Monte Carlo 模拟 结果 的 可 信和 度 较 
高 。 进 一 步 根据 线性 插值 法 求 得 统计 检验 力 为 
0.95 时 对 应 的 样本 量 约 为 488 人 。 


5 讨论 


51 ” 先 验 检验 力 分 析 时 的 参数 设 定 问题 
注意 到 ,人 研究 者 在 使 用 Monte Carlo 法 或 
Satorra-Saris 法 进行 先 验 检 验 力 分 析 时 , 均 需 要 指 
定理 论 协 方差 矩阵 。 该 步骤 中 ,不论 研究 者 是 采 
用 直接 设 定 协 方差 矩阵 的 方式 还 是 通过 设 定理 论 
模型 进而 生成 协 方差 矩阵 的 方式 ， 都 需要 设 定 大 
量 参数 ， 这 些 参数 的 设 定 会 对 结果 产生 一 定 影响 ， 
且 影 响 规律 较为 复杂 。 例 如 ， 有 研究 要 对 一 个 两 
维度 ， 每 个 维度 包含 4 个 题目 的 量 表 进 行 验证 性 
因子 分 析 。 假 设 为 两 维度 间 不 存在 相关 , H 
立 的 情况 下 两 维度 之 间 存 在 相关 系数 为 0.3 的 正 
相关 ， 进 行 先 验 检验 力 分 析 时 ， 当 研究 者 设 定 载 
荷 为 0.7.0.75.0.8 时 , 设 定 显著 性 水 平 为 0.05, 统 
计 检 验 力 为 0.95, 使 用 Satorra-Saris 法 计算 得 到 
的 计划 样本 量 分 别 为 527 人 、472 人 和 429 人 , 使 
用 Monte Carlo 法 计算 得 到 的 计划 样本 量 分 别 为 
521 人 、464 人 和 417 Ao XUW, 本 文 4.2 部 分 给 
出 的 例子 中 , 若 将 所 有 控制 变量 与 其 他 变量 的 相 
关系 数 从 0.2 变 为 0.1, 在 其 他 条 件 不 变 的 情况 下 ， 
使 用 Satorra-Saris 法 计算 得 到 的 计划 样本 量 为 352 
A, 使 用 Monte Carlo 法 计算 得 到 的 计划 样本 量 为 
233 人 。 而 当 分 析 对 象 是 中 介 效 应 或 有 调节 的 中 
介 效 应 这 类 由 参数 函数 (通常 是 乘积 ) 构 建 的 指标 
时 ,车 研究 者 使 用 理论 模型 生成 协 方差 矩阵 ， 即 
使 是 同样 大 小 的 指标 值 ， 不 同 的 参数 也 会 导致 分 
析 结 果 的 不 同 。 例 如 ， 对 简单 中 介 效 应 模型 而 言 ， 
设 定 自 变量 到 因 变 量 的 路 径 系数 为 0.2， 中 介 效 应 
值 为 0.16， 显 著 性 水 平 为 0.05， 统 计 检 验 力 为 
0.90， 中 介 前 半路 径 系 数 ( 自 变 量 到 中 介 变 量 的 路 
径 系数 ) 和 中 介 后 半路 径 系数 (中 介 变 量 到 因 变 量 
的 路 径 系 数 ) 分 别 取 0.2 和 0.8, 0.4 和 0.4, 0.8 F 
0.2( 所 有 系数 均 为 标准 化 系数 ) 的 情况 下 ， 使 用 
Satorra-Saris 法 计算 得 到 的 计划 样本 量 分 别 为 258 
人 、114 人 和 632 人 , 使 用 Monte Carlo 法 计算 得 
到 的 计划 样本 量 分 别 为 255 人 、88 人 和 620 人 。 


| 


chinaXiv:202303.09766v1 


ChinaXiv 合 作 期 刊 


2126 心理 科学 进展 第 30% 


从 上 述 例子 中 不 难看 出 先 验 检验 力 分 析 中 参数 设 
定 对 结果 的 影响 。 这 意味 着 , 研究 者 在 使 用 
Satorra-Saris 法 或 Monte Carlo 法 进行 先 验 检验 力 
分 析 时 ,应 当 和 审慎 确定 参数 的 大 小 ， 必 要 时 , 在 


错误 地 在 统计 检验 力 分 析 过 程 中 混入 了 抽样 误差 
而 降低 了 统计 检验 力 的 可 靠 性 ,并 进一步 提出 了 
纠正 该 误差 的 方案 (Liu & Wang, 2019)。 对 该 问题 
的 深入 讨论 已 超出 了 本 文 的 范畴 ， 有 兴趣 的 读者 


缺乏 既往 研究 或 理论 的 支持 下 , 人 研究 者 可 以 反复 
设 定 参数 的 值 ， 并 从 结果 中 选取 一 个 较 保守 的 样 
本 量 。 
52 ”理论 协 方差 矩阵 的 产生 方式 

如 前 所 述 ， 结 构 方 程 模 型 统计 检验 力 分 析 的 
Monte Carlo 法 和 Satorra-Saris 法 在 分 析 时 需要 指 
定理 论 协 方差 矩阵 ， 且 理论 协 方差 矩阵 的 指定 存 
在 两 种 不 同 模式 : 直接 指定 或 通过 指定 备 择 模型 
生成 。 虽然 从 分 析 过 程 和 分 析 结 果 上 二 者 可 能 3 
无 不 同 之 处 , 但 从 人 逻辑 上 讲 , 不 同 的 理论 协 方差 


可 以 进一步 阅读 Liu 和 Wang (2019) 的 研究 。 
5.3” 非 中 心 交 分 布 近似 带 来 的 误差 

Satorra-Saris 法 与 MacCallum 法 这 两 种 分 析 
方法 从 本 质 上 来 说 都 是 利用 非 中 心 x 分布 作为 当 
Hi 成 立时 样本 似 然 比 统计 量 Tu 分 布 的 近似 。 
Steiger 等 人 (1985) 指 出 ， 这 种 近似 当 且 仅 当 样本 
满足 似 然 函数 关于 分 布 的 基本 假设 、 样 本 量 不 大 
小 、 拟 合 函 数值 不 是 太 大 的 情况 下 成 立 。 诚 然 ， 
Satorra 等 人 (1991) 的 研究 进一步 指出 ， 当 样本 服 
从 正 态 分 布 且 理论 模型 与 实际 数据 的 偏离 程度 不 


矩阵 生成 方式 意味 着 对 结构 方程 建 模 过 程 中 误差 
来 源 的 不 同 理解 。Jak 等 人 (2021) 指 出 作为 衡量 理 
论 模 型 与 真实 数据 间 差 异 的 似 然 比 统计 量 Ty X 
际 上 包含 了 两 部 分 的 误差 : 随机 抽样 造成 的 抽样 
误差 和 理论 模型 与 真实 数据 不 匹配 的 误差 。 事实 
上 , 理论 模型 与 真实 数据 不 匹配 的 误差 还 可 以 再 
分 为 两 种 不 同 的 来 源 : 来 自 模型 内 部 参数 设 定 错 
误 以 及 模型 本 身 对 总 体 的 简化 。 当 研究 者 通过 指 
定 备 择 模 型 的 方式 生成 理论 协 方差 矩阵 时 ， 研 究 
者 实际 上 默认 了 模型 内 部 参数 设 定 错误 是 导致 误 
差 的 主要 来 源 ， 而 当 研究 者 通过 直接 指定 的 方式 
生成 理论 协 方差 矩阵 时 ， 则 将 这 两 种 不 同 来 源 的 
误差 “打包 ”在 了 一 起 。 事 实 上 , 经典 的 Satorra- 
Saris 法 与 MacCallum 法 之 间 也 具有 类 似 的 关系 ， 
经 典 的 Satorra-Saris 法 更 加 强调 来 自己 经 结构 化 
了 的 模型 内 部 的 误差 ，MacCallum 法 则 通过 模型 
拟 合 指数 将 模型 内 部 的 误差 和 模型 对 总 体 简 化 的 
误差 放 在 一 起 考察 。 这 种 细微 的 差别 意味 着 研究 
者 需要 在 面 对 不 同 检验 的 情况 下 灵活 选择 理论 协 
方差 矩阵 的 生成 方式 : 当 有 较 强 的 理论 或 既往 研 
究 指 向 误差 可 能 存在 于 模型 内 部 时 , 使 用 指定 备 
择 模 型 的 生成 方式 ; 而 当 研 究 者 对 这 些 信息 掌握 
程度 较 少 时 , 使 用 直接 指定 的 方式 或 干脆 换 用 
MacCallum 法 进行 分 析 。 

此 外 , 另 一 个 有 趣 的 议题 是 , 使 用 先前 的 研 
究 信 息 作 为 指定 理论 协 方差 的 依据 是 否 是 合适 
的 。 诚 然 ， 既往 人 研究 中 确实 存 有 使 用 相似 主题 的 
研究 结果 作为 统计 检验 力 分 析 过 程 中 设 定 参数 依 
止 的 先例 , 但 有 研究 者 指出 ， 这 种 指定 方式 可 能 


是 十 分 严重 时 ， 即 使 对 小 样本 而 言 ， 这 种 利用 非 
中 心 90 分 布 进行 近似 的 精度 也 是 基本 令 人 满意 
的 。 但 对 于 分 布 假 设 不 满足 时 这 种 近似 的 精度 是 
否 仍 令 人 满意 ， 目 前 仍 缺 乏 相关 研究 的 分 析 。 这 
意味 着 研究 者 在 对 结构 方程 模型 的 统计 检验 力 进 
行 分 析 时 ,应 当 充 分 了 解 样本 分 布 情 况 与 所 选 结 
构 方 程 模型 参数 估计 方法 之 间 的 适 配 性 ,在 进行 
统计 检验 力 分析 和 实际 进行 数据 处 理 时 使 用 同样 
的 估计 方法 , 在 样本 分 布 偏离 正 态 分 布 的 情况 下 ， 
也 可 以 考虑 换 用 适用 范围 更 广 的 Monte Carlo 法 
进行 统计 检验 力 分 析 。 目 前 主流 软件 提供 的 参数 
估计 方法 中 ,， 极 大 似 然 法 (maximum likelihood, 
ML) 和 广义 最 小 二 乘法 (generalized least squares, 
GLS) 要 求 数 据 为 正 态 分 布 , 加 权 最 小 二 乘法 
(weighted least squares, WLS)、 非 加 权 最 小 二 乘法 
(unweighted least squares, ULS) 和 对 角 线 加 权 最 小 
二 乘法 (diagonally weighted least squares, DWLS) 
对 数据 分 布 要 求 不 高 , 但 ULS 属于 尺度 依赖 方法 
(scale-dependent method), 使 用 时 最 好 保证 各 个 
变量 的 测量 单位 统一 , 而 WLS 和 DWLS 通 常 则 需 
要 1000 以 上 的 样本 量 ( 吴 明 隆 , 2010)。 

此 外 , Satorra-Saris 法 在 使 用 非 中 心 x 分布 对 
样本 似 然 比 统计 量 Tw 的 分 布 进行 近似 时 ， 同时 
使 用 模型 拟 合 “真实 总 体 ” 得 到 的 统计 量 作 为 对 
中 心 参数 的 估计 。 这 意味 着 使 用 Satorra-Saris X 
分 析 模 型 参数 检验 的 统计 检验 力 时 ， 检 验方 法 本 
身 也 会 对 结果 产生 影响 。 本 文 4.2 部 分 模型 参数 
检验 的 先 验 检验 力 分 析 中 ，Monte Carlo 法 与 
Satorra-Saris 法 对 同一 检验 的 分 析 得 出 了 差异 较 
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大 的 结果 就 是 这 个 原因 。 在 该 例 中 , 在 设 定 a = 
0.05, N = 401 的 情况 下 , 使 用 Wald 检验 、LM 检 
验 和 似 然 比 检验 得 到 的 统计 量 对 非 中 心 参数 进行 
估计 时 得 到 的 非 中 心 参数 分 别 为 13.399，22.658 
和 23.265， 对 应 的 统计 检验 力 分 别 为 0.955, 0.997 
和 0.998。 使 用 Monte Carlo 法 对 三 种 方法 进行 统 
计 检 验 力 分 析 ， 得 到 的 结果 则 分 别 为 0.996, 0.998 
和 0.998。 这 提示 研究 者 在 使 用 Satorra-Saris 法 分 
析 Wald 检验 的 统计 检验 力 时 ， 可 能 会 由 于 低估 统 
计 检 验 力 而 高 估 所 需要 的 样本 数量 。 
5.4 Monte Carlo 法 的 精度 与 计算 效率 

诚然 Monte Carlo 法 由 于 可 以 广泛 适用 于 各 
类 检验 场景 ， 且 在 模拟 次 数 足 够 多 的 情况 下 得 出 
的 统计 检验 力 的 估计 值 依 概率 收敛 于 其 理论 真 值 ， 
使 得 该 方法 广泛 受到 研究 者 的 青睐 。 但 Monte 
Carlo 法 同样 面临 着 结果 不 稳定 以 及 计算 效率 低 
下 的 情形 。 与 其 他 两 种 方法 不 同 , Monte Carlo 法 
得 到 的 计算 结果 并 不 是 一 个 确定 值 ， 而 是 一 个 随 
机 变量 ， 该 随机 变量 的 期 望 值 等 于 统计 检验 力 的 
真实 值 ， 这 就 导致 了 对 同一 个 假定 样本 量 进 行 多 
次 统计 检验 力 分 析 时 可 能 会 得 出 不 同 的 结果 。 这 
意味 着 人 研究 者 在 使 用 Monte Carlo 法 进行 统计 检 
验 力 分 析 时 求 得 的 样本 量 理应 为 一 约 数 而 非 确 切 
值 。 换 言 之 ,，Monte Carlo 法 无 法 像 Satorra-Saris 
法 或 MacCallum 法 那样 给 出 “理论 上 达到 目标 统 


重复 次 数 R = 1000, 样本 量 N = 401 WY, 使 用 
Monte Carlo 法 分 析 Wald 检验 统计 检验 力 所 需 的 
时 间 大 致 为 60 s (单线 程 ); 若 将 检验 方法 更 换 为 
LM 检验 或 似 然 比 检验 ， 执行 同样 的 分 析 所 需 时 
间 则 上 升 到 约 900 s (单线 程 ); A Xt AE BS BW 
Bootstrap 检验 的 统计 检验 力 进 行 分 析 ， 即 使 采用 
并 行 计 算 技 术 将 任务 分 配 到 多 线程 计算 机 的 不 同 
线程 上 处 理 ， 所 需 时 间 通 常 也 会 超过 一 个 小 时 。 
此 时 若 反复 设 定 样本 量 甚 至 参数 进行 多 轮 分 析 从 
时 间 上 来 说 是 不 经 济 的 。 诚 然 有 学 者 指出 此 时 可 
以 换 用 参数 Bootstrap 法 (或 与 之 等 价 的 MC 法 ) 简 
化 计算 过 程 ， 节 省 计算 时 间 (Schoemann et al., 
2017), 但 这 种 更 换 也 将 在 一 定 程度 上 造成 统计 检 
验 力 估计 时 的 误差 。 
5.5 不 同 拟 合 指数 的 选择 对 MacCallum 法 的 影响 
本 文 主要 介绍 了 用 MacCallum 法 对 RMSEA 
相关 的 零 假 设 进行 统计 检验 力 分 析 的 理论 与 实 
Eko 事实 上 , MacCallum 法 可 以 用 于 多 种 拟 合 指数 
的 情形 。MacCallum 和 Hong (1997) 将 该 方法 进 一 
步 扩展 到 了 针对 GFI 和 AGFI 相关 零 假 设 的 情形 。 
Kim (2005) 则 提出 了 当 零 假设 是 关于 RMSEA, 
CFI MFI 和 GFI 这 四 种 参数 其 一 时 估算 样本 量 的 
简便 公式 。 值 得 注意 的 是 ， 上 述 参数 间 的 换算 关 
系 还 会 受到 模型 本 身 的 影响 。 因 此 ， 目 前 学 界 主 
流 建议 的 拟 合 指数 界 值 在 不 同 的 模型 中 并 不 等 价 


NS 


计 检 验 力 所 需 的 最 小 样本 量 ”。 这 也 同时 意味 着 通 
过 缩短 采样 间距 的 方式 对 基于 Monte Carlo 法 的 
统计 检验 力 分 析 的 精度 的 提升 是 有 限 的 。 

此 外 , Monte Carlo 法 的 计算 效率 一 直 以 来 也 
是 比较 困扰 研究 者 的 一 个 问题 。 与 Satorra-Saris 
法 和 MacCallum 法 误差 不 同 , Monte Carlo 法 的 误 
差 来 源 主要 是 模拟 样本 带 来 的 偶然 误差 。 根据 De 
Moivre-Laplace 中 心 极限 定理 ， 随 着 重复 次 数 R 
的 提升 , Monte Carlo 法 所 得 到 的 统计 检验 力 估 计 
值 的 标准 差 会 降低 , 但 其 降低 速率 与 重复 次 数 R 
之 间 并 非 线 性 关系 ， 而 是 二 次 关系 。 换 言 之 ， 知 研 
究 者 希望 将 Monte Carlo 法 得 到 的 估计 值 的 精度 
提高 一 个 数量 级 (相当 于 让 估计 值 的 标准 差 降 低 
为 原来 的 0.1 倍 )， 则 需要 将 R 扩大 为 原来 的 100 
倍 。 这 大 大 降低 了 Monte Carlo 法 的 计算 效率 。 并 
H, Monte Carlo 法 本 身 计 算 起 来 就 比较 耗 时 , 在 
CPU 主 频 为 2.50 GHz 的 笔记 本 电脑 上 运行 本 文中 
关于 中 介 模 型 统计 检验 力 分 析 的 R 程序 ， 当 设 定 


这 就 导致 了 针对 同一 个 计划 模型 ， 采 用 不 同 拟 合 
指数 进行 先 验 检 验 力 分 析 时 ， 可 能 得 出 不 同 的 结 
果 。 例 如 ,根据 温 忠 刨 等 人 (2004) 的 建议 ， 评 佑 模 
型 拟 合 时 ，MFI 的 建议 界 值 为 0.8$，RMSEA 的 建 
议 界 值 为 0.08。 此 时 , 对 于 不 同 自 由 度 的 模型 而 
言 ， 依据 这 两 个 指标 的 建议 界 值 计算 出 的 计划 样 
本 量 也 有 所 不 同 。 以 计划 显著 性 水 平 为 0.05, 统 
计 检 验 力 为 0.95 的 情况 为 例 ， 当 模型 自由 度 为 10 
时 , 根据 MFI 和 RMSEA 计算 出 的 计划 样本 分 别 
有 77 人 和 383 A; 当 模 型 自由 度 为 50 时 ， 根 据 
MFI 和 RMSEA 计算 出 的 计划 样本 分 别 有 134 人 
和 136 A; 而 当 模 型 自由 度 为 100 时 , 根据 MFI 
和 RMSEA 计算 出 的 计划 样本 则 分 别 有 176 人 和 
90 人 。 因 此 , 本 文 建议 研究 者 在 使 用 MacCallum 
法 进行 先 验 检验 力 分 析 时 , 使 用 多 种 拟 合 指 数 构 
建 假设 进行 检验 ,并 在 所 有 计划 样本 量 中 选取 最 
大 值 。 至 于 拟 合 指数 的 具体 选择 ， 可 以 参考 温 忠 
麟 等 人 (2004) 的 研究 。 
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5.6 ”评价 模型 拟 合 时 的 统计 检验 力 分 析 与 等 效 

如 前 所 述 ， 既 往 研 究 中 对 传统 x 检验 的 一 个 
批评 就 是 该 检验 将 想 要 验证 的 假设 作为 零 假 设 ， 
在 结果 报告 中 又 缺乏 对 其 开 类 错误 发 生 概 率 的 分 
析 ， 使 其 无 法 真正 起 到 验证 模型 的 作用 ( 王 阳 等 ， 
2020)。 而 等 效 性 检验 实质 上 就 是 通过 构建 一 个 


性 水 平 为 B 的 水 平 下 ， 对 实际 拟 合 良 好 的 模型 ( 即 
Furo =), RAAF IBBON Ho: Furo > so 的 等 效 性 
检验 进行 检验 时 ,正确 拒 绝 零 假 设 的 概率 为 1-a。 
因此 ,本 文 建议 研究 者 在 对 评估 模型 拟 合 相关 的 
假设 检验 进行 统计 检验 力 分 析 时 ， 设 定 不 小 于 
0.95 的 统计 检验 力 。 这 样 可 以 保障 在 模型 拟 合 较 
差 的 情况 下 , 研究 者 错误 地 认为 模型 拟 合 良 好 的 


“可 接受 的 差异 ”标准 的 方式 重 写 零 假 设 ， 从 而 将 
想 要 验证 的 假设 重新 放 回 了 备 择 假 设 的 位 置 ， 在 
一 定 程度 上 ， 也 可 以 将 等 效 性 检验 视 作 一 种 “ 颠 
IT H 5 H, WIER x 检验 ”"。 根 据 MacCallum 
等 人 (1996) 的 研究 ， 可 以 推 知 , 对 Hoa: Furo S ko, 
且 真 实 参数 为 Furo = ky 的 情形 ， 存 在 一 种 与 之 
“对 偶 ?” 的 情形 : Ho: Furo > 局 ， 且 此 时 的 真实 参 
数 为 Pirro = Ko， 这 两 种 情形 下 I 类 错误 概率 和 开 
类 错误 概率 互相 为 对 方 的 荆 类 错误 概率 与 1 类 错 


可 能 性 不 高 于 0.05。 
6 总 结 


本 文 主要 介绍 了 结构 方程 模型 假设 检验 统计 
检验 力 分 析 的 三 种 代表 方法 : Satorra-Saris 法 、 
MacCallum 法 和 Monte Carlo 法 ， 这 些 方法 各 有 自 
己 的 适用 情形 。 相 比 于 目前 已 被 广大 研究 者 所 熟 
知 的 平均 数 差 异 检验 的 统计 检验 力 分 析 ，, 结构 方 
程 模型 的 统计 检验 力 分 析 对 研究 者 的 统计 水 平 有 


误 概 率 。 这 意味 着 ， 对 传统 x 检验 而 言 ， 在 进行 
先 验 检验 力 分 析 时 ， 若 研究 者 参考 等 效 性 检验 中 
的 Hot Bee Fur 的 真 值 ( 即 Fy 的 真 值 不 小 于 so)， 并 

设 定 显著 性 水 平 为 a, 统计 检验 力 为 1-B 计算 出 
计划 样本 量 N。 那 么 N 同样 可 以 保证 在 控制 显著 


基于 x? 分布 的 方法 基于 模拟 或 重 抽样 的 方法 Æ 


aa 假设 检验 的 目的 是 ?| | Monte Carlo 法 


更 高 的 要 求 ， 并且 检验 过 程 较为 复杂 。 综 合 考 虑 
易 掌握 和 易 使 用 两 方面 的 要 求 ， 本 文 建议 研究 者 
在 研究 设计 的 过 程 中 根据 自己 需要 进行 的 检验 类 
型 选择 恰当 的 统计 检验 力 分 析 方式 ， 具体 可 以 参 
考 图 2 进行 选择 。 


先 验 检验 力 分 析 还 是 
后 验 检验 力 分 析 ? 


计划 使 用 的 是 否 为 基于 ?分 布 
假设 检验 方法 是 ”| | 的 模型 拟 合 检验 ? 


i 
检验 对 PT 
是 否 有 明确 备 择 模型 ”| | (如 路 径 系数 ) 还 是 复 
een 中 介 效 应 ) 
REM 
mwen 
ETA 计划 采用 的 假设 
理论 协 方差 矩阵 的 = Satorra-Saris 法 检验 方法 是 
Satorra-Saris 法 r LM 检验 或 似 然 比 检验 
Wald 检 验 
y 
Monte Carlo 法 Satorra-Saris 法 


图 2 ”结构 方程 模型 统计 检验 力 分 析 方 法 选择 建议 
注 : 未 标注 理论 协 方差 矩阵 生成 方法 的 表示 研究 者 可 以 参照 5.2 部 分 的 建议 自行 设 定 
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Power analysis in structural equation modeling: Principles and methods 


ZHAI Hongkun, LI Qiang, WEI Xiaowei 
(Department of Social Psychology, Nankai University, Tianjin 300350, China) 


Abstract: Structural equation modeling (SEM) is an important statistical tool in psychology, management, 
and sociology. However, the lack of analyses and reports of statistical power in many studies using SEM has 
reduced the probative power of the study results. There are three main types of statistical power analysis 
methods for SEM: the Satorra-Saris, MacCallum, and Monte Carlo methods. The Satorra-Saris method 
applies to cases where the alternative model is clear, analysis object is simple, and the statistic to be 
examined is based on the x° distribution. The MacCallum method is applicable to the case of y°-based 
overall model fit tests with unknown alternative models. Further, the Monte Carlo method applies to cases 
where the analysis object is complex, or the test is performed using simulation or resampling methods. In 
practice, researchers should clarify the purposes and methods of the test as well as the availability of 
alternative models as a priority. The research method can then be determined on the basis of the above 
information. 


Key words: structural equation model, statistical power, fitting test, test of model parameters 
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《结构 方程 模型 统计 检验 力 分 析 : 原理 与 方法 》 补 充 材 料 


Ostet A: 案例 部 分 的 R 代码 

# 载 人 程序 包 

if(!require(lavaan) ){install.packages('lavaan')} 
if(!require(mvtnorm) ) {install.packages('mvtnorm' )} 
if (!require(semPower )){install.packages('semPower' )} 


#4.1 模型 拟 合 检验 的 先 验 检验 力 分 析 
#4.1.1 Satorra-Saris 法 
#(1) 设 定 用 于 生成 理论 协 方差 矩阵 的 模型 


genMod = ' 

POS ~~ 0.23 * DES 

POS ~~ 0.40 * ANG 

DES ~~ 0.78 * ANG 

POS =~ 0.87 * X1 + 0.90 * X2 + 0.75 * X3 + 0.55 * X4 
DES =~ 0.59 * X5 + 0.72 * X6 + 0.58 * X7 + 0.71 * X8 + 0.59 * X9 
ANG =~ 0.73 * X10 + 0.78 * X11 + 0.57 * X12 

X1 ~~ 0.24 * X1 

X2 ~~ 0.19 * X2 

X3 ~~ 0.44 * X3 

X4 ~~ 0.70 * X4 

X5 ~~ 0.65 * X5 

X6 ~~ 0.47 * X6 

X7 ~~ 0.54 * X7 

X8 ~~ 0.49 * X8 

X9 ~~ 0.52 * X9 


X10 ~~ 0.45 * X10 

X11 ~~ 0.39 * X11 

X12 ~~ 0.55 * X12! 

#(2) 设 定 用 于 分 析 的 模型 

anaMod = ' 

POS =~ X1 + X2 + X3 + X4 

DES =~ X5 + X6 + X7 + X8 + X9 
ANG =~ X10 + X11 + X12' 

#(3) 生成 理论 协 方差 矩阵 

E1 = cfa(genMod)@implied$cov[ [1] ] 
rownames(E1) = pasteO('X', 1:12) 
#(4) 分 析 模型 

fit = cfa(anaMod, sample.cov = E1, sample.nobs = 101, likelihood = 'Wishart' ) 
#(5) 提取 关键 参数 

fo = fit@optim$fx * 2 

#(6) 用 于 计算 等 效 性 检验 样本 量 的 函数 
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setPowerFindEqTestN = function(f0, df, alpha = 0.05, beta = 0.8, e = 0.08){ 
foo = function(N){pchisq(qchisq(alpha, df, ncp = e*e*df*N), df, fO*N) - beta} 
n = uniroot(foo, c(10, 1e4))$root 
return(ceiling(n+1) ) 

} 

#(7) 执行 计算 

setPowerFindEqTestN(f0, df = 51, alpha = 0.05, beta = 0.95) 


#4.1.2 MacCallum 7: 

semPower.aPriori(effect = 0.08, effect.measure = 'RMSEA', alpha = 0.05, 
power = 0.95, df = 51) -> RES 

summary (RES) 


#4.1.3 Monte Carlo 法 
#(1) 生成 理论 协 方差 矩阵 的 部 分 与 Satorra-Saris HA, E HIEN TT 
#(2) 用 于 计算 统计 检验 力 的 程序 
mcCFA = function(N, E, mod, R = 1000L){ 
RES = matrix(0, nrow = R, ncol = 62) 
for(i in 1:R){ 
Sample = rmvnorm(N, sigma = E) 
if(!is.null(colnames(E))){ 
colnames(Sample) = colnames(E) 
selse{ 
colnames(Sample) = rownames(E) 


} 
RES[i, ] = tryCatch({ 
tempFit = cfa(mod, Sample, std.lv = TRUE, likelihood = 'wishart') 
ind = sum(unlist(parTable(tempFit)[, 14:15])) 
ind = ind / ind 
c(!is.nan(ind), unlist(parTable(tempFit)[, 14:15]), 
as.numeric(fitmeasures(tempFit, 'chisq'))) 


} 

, warning = function(w){c(0, rep(NaN, 61))}) 
} 
return(RES) 


} 

#(3) 利用 程序 计算 统计 检验 力 

set.seed(1234) 

mcCFA(101, E1, anaMod) -> RES1 

mcCFA(201, E1, anaMod) -> RES2 

#(4) 统计 结果 ， 并 进行 规范 化 输出 

resCFA = function(RES, TV, N){ 
R = nrow(RES) 
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2B 


P = sum(RES[, 1]) 7 R 

RES[which(RES[, 1] == 1), -1] 

COV = matrix(@, nrow = nrow(RES), ncol = 30) 

for(i in 1:nrow(RES)){for(j in 1:30){COV[i, j] = (TV[j] <= qnorm(0.975, RES[i, j], 
RES[i, j+30])) & (TV[j] >= qnorm(®0.025, RES[i, j], RES[i, j+30]))}} 

ME1 = apply(RES[, 1:30], 2, mean) 

ME2 apply(RES[, 31:60], 2, mean) 

SE = apply(RES[, 1:30], 2, sd) 

PC = apply(COV, 2, mean) 

PM = cbind(ME1, ME2, SE, PC, abs(ME1 - TV) / TV, abs(ME2 - SE) / SE) 

RES = RES[, 61] 

Power = sum(RES <= qchisq(0.05, 51, 0.08*0.08*51*(N-1))) / length(RES) 

colnames(PM) = c('MC_MEAN', 'MC_SE', 'MEAN_SE', 'Coverage', 'Mean_Bias', 

"SE _Bias') 
RES = list(Proper = P, CHECK = PM, POWER = Power) 
return(RES) 


Pas] 
m 
n 
ll 


} 

#(5) 给 出 参数 的 真实 值 

TV = parTable(cfa(anaMod, sample.cov = E1, sample.nobs = 101, 
std.lv = TRUE, likelihood = 'Wishart'))[, 14] 

#(6) 给 出 规范 化 输出 

resCFA(RES1, TV, 101) 

resCFA(RES2, TV, 201) 

#(7) 补充 中 点 采样 

mcCFA(151, E1, anaMod) -> RES3 

resCFA(RES3, TV, 151) 


#42 模型 参数 检验 的 先 验 检验 力 分 析 
#4.2.1 Satorra-Saris 法 ; Wald 检验 
#(1) 设 定理 论 协 方差 矩阵 


E2 = matrix(c( 


1.0, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 
0.1, 1.0, 0.1, 0.1, 0.1, 0.1, 0.1, 
0.1, 0.1, 1.0, 0.1, 0.1, 0.1, 0.1, 
0.1, 0.1, 0.1, 1.0, 0.3, 0.3, 0.3, 
0.1, 0.1, 0.1, 0.3, 1.0, 0.3, 0.3, 
0.1, 0.1, 0.1, 0.3, 0.3, 1.0, 0.3, 
0.1, 0.1, 0.1, 0.3, 0.3, 0.3, 1.0 


), nrow = 7) 

colnames(E2) = c('Age', 'Ser', 'STR', 'ECI', 'PSS', 'RES', 'SRG') 
#(2) 设 定 分 析 模 型 

anaMod = ' 

SRG ~ b*RES + d*PSS + ECI + STR + Ser + Age 
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RES ~ k*PSS + ECI + STR + Ser + Age 
PSS ~ a*ECI + e * STR + Ser + Age 
akb := a*k*b' 
#(3) 用 于 计算 统计 检验 力 的 程序 
SSSEM = function(E, mod, PID, alpha = 0.05, beta = 0.95){ 
foo = function(N){ 
fit = sem(mod, sample.cov = E, sample.nobs = N, likelihood = 'Wishart') 
NCP = (parTable(fit)[PID, 14] / parTable(fit)[PID, 15]) ** 2 
return(pchisq(qchisq(1 - alpha, 1), 1, NCP, FALSE) - beta) 
} 
n = uniroot(foo, c(10, 1e4))$root 
return(ceiling(n + 1)) 
} 
#(4) 执行 计算 
ssSEM(E2, anaMod, 29) 


#4.2.2 Monte Carlo 法 ; wald 检验 
#(1) 生成 理论 协 方差 矩阵 的 部 分 与 Satorra-Saris 法 相同 ， 这 里 直接 使 用 了 前 面 的 结果 
#(2) 用 于 计算 统计 检验 力 的 程序 
mcSEM = function(N, E, mod, R = 1000L){ 
# 注 : 基于 无 完全 共 线 性 的 真实 数据 协 方差 矩阵 的 显 变量 回归 理论 上 不 可 能 出 现 不 恰当 解 ， 因 此 没有 考虑 对 其 
进行 统计 
RES = matrix(0, nrow = R, ncol = 58) 
for(i in 1:R){ 
Sample = rmvnorm(N, sigma = E) 
if(!is.null(colnames(E))){ 
colnames(Sample) = colnames(E) 
selse{ 


colnames(Sample) = rownames(E) 


} 

for(j in 1:ncol(Sample)){Sample[, j] = (Sample[, j] - mean(Sample[, j])) / 
sd(Sample[, j])} 

tempFit <- sem(mod, Sample, likelihood = 'wishart') 

RES[i, ] = unlist(parTable(tempFit)[, 14:15]) 


I 
return(RES) 


} 

#(3) 利用 程序 计算 统计 检验 力 
set.seed(4321) 

mcSEM(601, E2, anaMod) -> RES1 
mcSEM(501, E2, anaMod) -> RES2 
mcSEM(401, E2, anaMod) -> RES3 
mcSEM(301, E2, anaMod) -> RES4 
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#(4) 统计 结果 ， 并 进行 规范 化 输出 
resSEM = function(RES, TV){ 
R = nrow(RES) 
COV = matrix(@, nrow = nrow(RES), ncol = 29) 
for(i in 1:nrow(RES)){for(j in 1:29){COV[i, j] = (TV[j] <= qnorm(0.975, RES[i, j], 
RES[i, j+29])) & (TV[j] >= qnorm(0.025, RES[i, j], RES[i, j+29]))}} 
ME1 = apply(RES[, 1:29], 2, mean) 
ME2 = apply(RES[, 30:58], 2, mean) 
SE = apply(RES[, 1:29], 2, sd) 
PC = apply(COV, 2, mean) 
PM = cbind(ME1, ME2, SE, PC, abs(ME1 - TV) / TV, abs(ME2 - SE) / SE) 
Power = sum(abs(RES[, 29] / RES[, 58]) >= qnorm(0.975)) / R 
colnames(PM) = c('MC_MEAN', 'MC_SE', 'MEAN_SE', 'Coverage', 'Mean_Bias', 
"SE _Bias') 
RES = list(CHECK = PM[1:15, ], POWER = Power) 
return(RES) 


} 

#(5) 给 出 参数 的 真实 值 

TV = parTable(cfa(anaMod, sample.cov = E2, sample.nobs = 101, 
likelihood = 'Wishart'))[, 14] 

#(6) 给 出 规范 化 输出 

resSEM(RES1, TV) 

resSEM(RES2, TV) 

resSEM(RES3, TV) 

resSEM(RES4, TV) 

# 结果 显示 N = 301/601 时 精度 不 足 

#(7) 重新 抽取 精度 不 足 的 样本 

mcSEM(601, E2, anaMod) -> RES1 

mcSEM(301, E2, anaMod) -> RES4 

resSEM(RES1, TV) 

resSEM(RES4, TV) 


@ 补 充 材 料 B: 统计 检验 力 的 后 验 分 析 实 例 

本 部 分 的 演示 样 例 同样 来 自 Zhai 等 人 (2021) 的 研究 。 该 研究 样本 量 为 439 人 ,其 中 男 157 A, 女 282 
人 。, 本 部 分 同样 以 该 研究 中 对 情绪 调节 自我 效能 感 量 表 的 验证 性 因子 分 析 与 以 压力 后 成 长 为 因 变量 的 链 
式 中 介 模 型 为 分 析 对 象 ， 对 其 进行 后 验 统计 检验 力 分 析 。 该 研究 中 ,关于 情绪 调节 自我 效能 感 量 表 的 验 
证 性 因子 分 析 结 果 为 : y = 151.979, df= 51, RMSEA = 0.067, 参考 王 阳 等 人 (2020) 的 建议 ,以 RMSEA, = 
0.08 为 经 验 临界 值 进行 等 效 性 检验 得 p = 0.046。 链 式 中 介 模 型 所 涉及 变量 间 的 协 方差 矩阵 参见 表 B-1, 
具体 的 模型 设置 同 正文 部 分 ， 路径“ 情绪 创造 力 一 领悟 社会 支持 一 情绪 调节 自我 效能 感 一 压力 后 成 长 ” 
中 介 效 应 值 为 0.009( 标 准 化 结果 为 0.021), Wald 检验 z= 2.713, p= 0.007。 参 考 正文 图 2 中 的 建议 , 综合 
确定 使 用 MacCallum 法 对 模型 拟 合 检验 部 分 进行 后 验 检验 力 分 析 , 使 用 Monte Carlo 法 对 链 式 中 介 效 应 
部 分 进行 后 验 检验 力 分 析 。 


HI 
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表 B-1. 链 式 中 介 模 型 所 涉及 变量 间 的 协 方差 矩阵 
1 2 3 4 5 6 7 
1. 年 龄 36.818 
2. 疫情 严重 程度 —0.194 0.936 
3. 疫情 相关 压力 -0.643 2.788 215.412 
4. 情绪 创造 力 一 9.107 0.098 9.321 129.108 
5. 领悟 社会 支持 -3.390 一 0.992 0.305 60.613 100.388 
6. 情绪 调节 自我 效能 感 一 1.217 一 0.294 2.198 42.998 29.468 31.922 
7. 压力 后 成 长 一 1.981 -0.275 —6.180 


; é $ i 24.855 21.677 11.936 21.294 


1. 模型 拟 合 检验 的 后 验 检 验 力 分 析 


在 后 验 检验 力 分 析 中 , 通常 直接 以 得 出 的 结果 作为 真 值 。 故 本 例 相 当 于 检验 村 
RMSEA 真 值 为 0.067 时 , 样本 的 入 统计 量 落 在 等 效 性 检验 零 假 设 Hy: RMSEA > 


Th 


本 量 为 439、a = 0.05, 
0.08 拒绝 域 中 的 概 


率 。 根 据 正文 中 的 公式 (1a) 和 (2a)， 可 以 求 出 后 验 检 验 力 ， 对 应 的 R 脚本 如 下 : 
= (439-1) * 51 * 0.080**2 
= (439-1) * 51 * 0.067**2 
= qchisq(0.05, 51, ncpO) 
= pchisq(Tmlc, 51, ncpa)) 
— 可 得 后 验 检验 力 为 0.546。 
.中 介 效 应 检验 的 后 验 检验 力 分 析 


ncpo 
ncpa 
Tmlc 
(POW 


ae ennai 
检验 为 分 析 对 象 ,控制 1 类 错 


中 关键 参数 的 估计 误差 均 不 超过 5%, 


o 这 主要 是 由 于 这 些 参数 的 理论 
这 些 参数 都 是 与 控制 变 
ee 


RERE 


"情绪 调节 上 
BIRKE a 为 0.05， 


我 效能 感 一 压力 后 成 长 "这 一 路 径 中 介 效 应 的 Wald 
WIE N 为 439 进行 分 析 ， 得 到 的 结果 全 部 为 恰当 解 ， 对 
应 的 统计 检验 力 为 0.847， 对 应 的 各 参数 佑 计 情 况 、 相 对 误差 以 及 95% 和 置信 区 间 覆 盖 率 参见 表 B-2。 本 例 
95% 置 信 区 间 覆 盖 率 均 处 于 0.91 至 0.98 之 间 ， 部 分 参数 的 均值 误 
值 过 小 导致 的 ， 


实际 上 这 些 参 数 标 准 误 的 误差 均 不 超过 5%, 


量 相关 的 参数 ， 说 明 Monte Carlo 模拟 结果 仍 具 可 信和 度 。 这 部 分 的 R 脚本 为 : 


if (!require(lavaan) ) {install.packages('lavaan' )} 
# 读 和 样本 协 方差 矩阵 


E = matrix(c( 


29,390, -1.217, -1.981, 
0.098, -0.992, -0.294, -0.275, 
0.305, 2.198, -6.180, 
60.613, 42.998, 24.855, 
60.613, 100.388, 29.468, 21.677, 


42.998, 29.468, 31.922, 11.936, 
21.677, 11.936, 21.294), nrow = 7) 


36.818, -0.194, -0.643, -9.107, 

-0.194, 0.936, 2.788, 

-0.643, 2.788, 215.412, 9.321, 

-9.107, 0.098, 9.321, 129.108, 

-3.390, -0.992, 0.305, 

-1.217, -0.294, 2.198, 

-1.981, -0.275, -6.180, 24.855, 
rownames(E) = c('Age', 'Ser', 'STR', 'ECI', 
# 设 定 分 析 模型 

anaMod = ' 


SRG ~ b*RES + d*PSS + ECI + STR + Ser + Age 


RES ~ k*PSS + ECI + STR + Ser + Age 
PSS ~ a*ECI + e * STR + Ser + Age 


'PSS', 'RES', 'SRG') 
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akb := a*k*b' 
TV = parTable(sem(anaMod, sample.cov = E, sample.nobs = 439, likelihood = 'Wishart'))[, 14] 
# 执 行 分 析 ， 这 里 利用 wishart 分 布 奉 代 了 原始 数据 


mcSEM = 
RES 
SaS 


function(N, E, mod, R = 1000L){ 
matrix(0, nrow = R, ncol = 58) 
= rWishart(R, N-1, E) / (N-1) 


for(i in 1:R){ 


tempE = SaS[, , i] 

rownames(tempE) = rownames(E) 

tempFit = sem(mod, sample .nobs =N, sample .cov = tempE, likelihood = 'wishart') 
RES[i, ] = unlist(parTable(tempFit)[, 14:15])} 


return(RES)} 
set .seed(1234) 
RES = mcSEM(439, E, anaMod) 
resSEM = function(RES, TV){ 


R= 


nrow(RES) 


COV = matrix(0, nrow = nrow(RES), ncol = 29) 
for(i in 1:nrow(RES)){for(j in 1:29){ 


cov[i, j] = (TV[j] <= gnorm(0.975, RES[i, j], RES[i, j+29])) & 
(TV[j] >= qnorm(0.025, RES[i, j], RES[i, j+29]))}} 


ME1 = apply(RES[, 1:29], 2, mean) 
ME2 = apply(RES[, 30:58], 2, mean) 


SE 
PC 
PM 


apply(RES[, 1:29], 2, sd) 
apply(COV, 2, mean) 
cbind(ME1, ME2, SE, PC, abs(ME1 - TV) / TV, abs(ME2 - SE) / SE) 


Power = sum(abs(RES[, 29] / RES[, 58]) >= qnorm(0.975)) / R 
colnames(PM) = c('MC_MEAN', 'MC_SE', 'MEAN_SE', 'Coverage', 'Mean_Bias', 


"SE Bias') 


RES = list(CHECK = PM[1:15, ], POWER = Power) 
return(RES) } 
resSEM(RES, TV) 


表 B-2. 中 介 效 应 Monte Carlo 检验 的 参数 估计 情况 


路 径 理论 值 估计 均值 ”均值 标准 误 标准 误 均 值 覆盖 率 均值 误差 ”标准 误 误 差 
RES—SRG 0.145 0.144 0.045 0.045 0.949 0.002 0.002 
PSS 一 SRG 0.118 0.118 0.022 0.022 0.953 0.002 0.023 
ECI 一 SRG 0.090 0.090 0.023 0.023 0.945 0.003 0.006 
STR 一 SRG —0.034 —0.033 0.013 0.012 0.958 —0.012 0.013 
Ser 一 SRG —0.035 一 0.047 0.193 0.195 0.949 —0.324 0.009 
Age 一 SRG —0.017 —0.018 0.030 0.031 0.945 —0.065 0.023 
PSS 一 RES 0.126 0.127 0.023 0.022 0.961 0.009 0.035 
ECI 一 RES 0.277 0.278 0.020 0.021 0.941 0.002 0.014 


STR 一 RES 0.001 0.001 0.013 0.014 0.943 0.574 0.010 
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路 径 理论 值 。 ”估计 均值 。 均值 标准 误 ERS BHAR 。 ”均值 误差 ”标准 误 误 差 
Ser—>RES —0.203 一 0.201 0.205 0.201 0.950 —0.007 0.018 
Age—RES 0.046 0.046 0.032 0.032 0.949 0.007 0.011 
ECI~PSS 0.472 0.470 0.036 0.036 0.944 0.004 0.008 
STR 一 PSS —0.005 —0.004 0.028 0.029 0.942 —0.126 0.021 
Ser—PSS 一 1.091 一 1.091 0.423 0.431 0.941 0.000 0.017 
Age 一 PSS 0.019 0.023 0.067 0.068 0.946 0.198 0.021 


TE: 加 粗 的 表示 该 路 径 系数 是 关键 参数 ; RES: 情绪 调节 自我 效能 感 , PSS: 领悟 社会 支持 , ECI: 情绪 创造 力 , SRG: HE 
力 后 成 长 , STR: 疫情 相关 压力 , Ser: 疫情 严重 程度 , Age: 年 龄 
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@ 补 充 材 料 C: MacCallum 法 适用 的 拟 合 指数 及 其 与 非 中 心 参数 的 换算 公式 
根据 正文 部 分 的 论述 , 不 难 发 现 MacCallum 法 的 核心 在 于 将 拟 合 指数 根据 其 定义 式 逆 变换 为 非 中 
心 参数 ， 作 为 对 6 的 估计 值 。 同 时 , 根据 正文 中 的 (3) 式 , 不 难 有 : 


6=(N-1)Fno=(N Dm L) Ty, -df (C.1) 


换言之 ,， 如 果 研 究 者 可 以 将 拟 合 指数 “还 原 ” 为 原始 的 9? 值 ， 并 且 知 晓 模型 自由 度 ， 那么 理论 上 研究 
者 就 可 以 根据 该 参数 求 出 非 中 心 参数 6 的 估计 值 。 侯 杰 泰 等 人 (2004) 的 书 中 给 出 了 绝 大 多 数 结构 方程 模 
型 拟 合 指数 的 定义 ， 从 中 不 难看 出 ， 当 知晓 信息 足够 多 的 情况 下 ,除了 少数 直接 基于 模型 隐 含 协 方差 矩 
阵 定义 的 拟 合 指数 (如 RMR、SRMR 等 ) 外 ,其 余 的 拟 合 指数 均 可 以 “还 原 ” 为 原始 的 x? 值 并 进一步 计算 出 
非 中 心 参数 8。 但 不 同 拟 合 指数 的 使 用 频率 不 同 , 将 其 换算 为 非 中 心 参数 所 需要 的 额外 参数 (或 者 说 , 换 
算 的 “麻烦 程度 ”) 也 不 同 ,， 故 并 不 是 所 有 的 拟 合 指 数 都 在 实践 中 被 应 用 ,目前 较为 常用 的 拟 合 指 数 有 : 
RMSEA, MFI, GFI, AGFI 和 CFI。 这 里 特别 需要 注意 的 是 , 目前 软件 和 论文 中 对 GFI 和 AGFI 的 缩写 
存在 一 定 混用 的 情况 ， 可 能 不 加 区 别 地 将 “Goodness of fit index/LISREL GFI”( 下 称 GFID 和 “Unbiased 
goodness of fit index/Steiger’s y” (下 称 GFI”); “Adjusted goodness of fit index”( 下 称 AGFID 和 “Unbiased 
adjusted goodness of fit index”( 下 称 AGFI WEA. WAR 为 例 , R 语言 结构 方程 建 模 程序 包 lavaan 中 关于 
GFI 的 定义 采用 了 GFI, 而 semPower 包 中 则 采用 了 GFT 。 特 别 需 要 说 明 的 是 ，MacCallum 法 中 使 用 的 
是 GFT 和 AGEFT。 因 此 若 对 lavaan 中 得 到 的 GFI 值 进行 后 验 检 验 力 分 析 ， 则 可 能 出 现 不 正确 的 结果 。 这 
五 种 指数 与 非 中 心 参数 换算 的 公式 为 : 


D 


5 =(N-1)-df RMSEA? (C.1a) 

5 =(N-1)-(-2InMFI) (C.1b) 
p(1—GET’) 

5=(N-1) o (C.1c) 


pa- AGET )df 
p(p +1) —2df (1— AGFT ) 
ô =(1—CFD[(N -DF yun 一 dy] (C.le) 


ô=(N 


(C.1d) 
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其 中 ,p 为 协 方差 矩阵 的 秩 ( 即 变量 维度 ), Fww 和 df 分 别 为 虚 模 型 (Null Model) 的 拟 合 函数 值 和 自由 度 ， 
FF 中 Fy = —In|R p R 是 变量 的 相关 和 矩阵 。 
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补充 材料 D: 基于 G*Power 的 结构 方程 模型 统计 检验 力 分 析 

本 文正 文部 分 在 介绍 结构 方程 模型 统计 检验 力 分 析 时 主要 使 用 了 R 语言 作为 软件 工具 。 但 R 语言 
本 身 学 习 成 本 较 高 ,对 于 不 熟悉 编程 的 研究 者 而 言 不 是 非常 友好 。 事 实 上 , R 语言 并 非 进行 结构 方程 模 
型 统计 检验 力 分 析 的 唯一 途径 ， 对 于 那些 不 熟悉 R 的 研究 者 ， 只 要 研究 者 本 人 能 较为 熟练 地 使 用 目前 市 
面 上 常见 的 结构 方程 建 模 软件 (如 LISREL, EQS, AMOS, Mplus 等 )， 再 结合 目前 学 界 常用 的 G*Power, 
即 可 完成 基于 x 分 布 的 结构 方程 模型 统计 检验 力 的 分 析 ， 这 在 一 定 程度 上 也 能 够 帮助 研究 者 应 对 一 些 
不 是 那么 复杂 的 应 用 场景 。 

1. 使 用 G*Power 对 结构 方程 模型 统计 检验 力 进行 分 析 的 原理 

正如 正文 中 提 到 的 那样 ， 基 于 x 分 布 的 统计 检验 力 分 析 的 核心 包括 两 部 分 : 四 根据 一 些 方式 估计 出 
非 中 心 参数 ; @ 根 据 非 中 心 参 数 结合 公式 (19) 或 (1b) 算 出 对 应 情况 下 的 统计 检验 力 。 其 中 四 的 估计 又 可 分 
为 基于 协 方差 矩阵 理论 真 值 的 Satorra-Saris 法 和 基于 拟 合 指数 的 MacCallum 法 ， 这 部 分 是 需要 研究 者 自 
己 套 用 公式 计算 或 借助 结构 方程 建 模 软 件 进行 计算 的 部 分 。G*Power 主要 针对 的 是 步骤 @ 中 的 内 容 ， 即 
在 已 知 零 假设 成 立时 对 应 的 非 中 心 参数 5。、 备 择 假设 成 立时 对 应 的 非 中 心 参数 61、 自 由 度 df、 工 类 错 
误 水 平 a 时 ,根据 公式 (2a) 或 (2b) 先 计算 出 拒绝 域 的 临界 值 ， 再 根据 公式 (1a) 或 (1b) 计 算出 统计 检验 力 的 
过 程 。 

2. G*Power 中 的 Generic x test 

不 难看 出 步骤 @ 的 核心 是 计算 x 分 布 的 分 位 数 和 单 尾 累积 概率 , 在 G*Power 中 ,相关 功能 可 以 在 
Generic X2 test 中 找到 。 打 开 G*Power 后 ， 先 选择 检验 族 (Test Family) y? tests， 而 后 便 可 在 检验 方式 
(Statistical test) 一 栏 中 找到 Generic y? test, G*Power 的 Generic x test 部 分 主要 提供 了 四 个 工 
具 :Compromise, Criterion, Post hoc 和 Sensitivity。 其 中 后 三 个 工具 与 结构 方程 模型 统计 检验 力 分 析 有 关 ， 
它们 的 作用 分 别 是 : Criterion: 给 出 统计 检验 力 1-B、 非 中 心 参数 CE G*Power 中 用 入 表示 ) 和 自由 度 df, 
求 出 对 应 的 工 类 错误 水 平 a; Post hoc: 给 出 非 中 心 参 数 5、 工 类 错误 水 平 a 和 自由 度 df 求 出 对 应 的 统 
计 检 验 力 1-B; Sensitivity: 给 出 I 类 错误 水 平 a、 统 计 检 验 力 1-B 和 自由 度 df, 求 出 对 应 的 非 中 心 参数 。 
但 特别 需要 注意 的 是 ，G*Power 提供 的 这 些 工 具 都 默认 零 假 设 成 立时 对 应 的 非 中 心 参 数 5。= 0， 即 零 假 
设 成 立时 统计 量 的 分 布 为 中 心 化 x 分 布 ,在 结构 方程 模型 假设 检验 的 经 典 方法 中 ,除去 等 效 性 检验 外 的 
其 他 方法 均 能 满足 G*Power 的 默认 条 件 ， 因 此 可 以 直接 借用 G*Power 进行 分 析 。 等 效 性 检验 的 情况 较 
为 复杂 ， 背 等 效 性 检验 的 6, =0 时, 可 以 借用 正文 5.6 部 分 的 结论 , 构建 与 其 对 偶 情形 ， JF ELK u 与 B m 
倒 代 入 求解 ; 若 等 效 性 检验 的 51 > 0 时， 则 需 借 用 中 心 化 分布 为 “ 桥 ”， 建立 起 两 个 非 中 心 Ja hY 
联系 。 

3. 具体 操作 方法 与 示例 

如 前 所 述 ,， 基于 x? 分 布 的 统计 检验 力 分 析 的 第 一 部 分 主要 需要 研究 者 自己 套用 公式 计算 或 借助 结 
构 方 程 建 模 软件 进行 计算 ,本 文 对 此 步 又 不 进行 详 述 , 仅 简单 介绍 本 步 中 需要 注意 的 几 个 要 点 。 中 本 步 
中 的 难点 在 于 选用 Satorra—Saris 法 时 协 方差 矩阵 理论 真 值 的 生成 ,特别 是 基于 模型 的 理论 协 方差 矩阵 真 
值 的 生成 ， 这 需要 研究 者 对 结构 方程 模型 的 基本 原理 有 较 好 的 理解 水 平 。 目 前 市 面 上 大 部 分 结构 方程 建 
模 的 教材 都 对 这 部 分 进行 了 相应 的 介绍 ， 系 统 学 习 过 结构 方程 建 模 的 研究 者 通常 可 以 在 教材 指导 下 自 
己 完 成 理论 协 方差 矩阵 的 生成 工作 。 另 外 , 生成 理论 协 方差 矩阵 时 需要 计算 和 矩阵 乘法 ， 这 部 分 可 通过 
Excel 的 MMULT 函数 完成 。@) 不 同 软件 计算 邓 统 计量 的 方法 略 有 不 同 ， 对 Mplus 而 言 , x = NFyr， 而 对 
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LISREL EQS 和 AMOS 而 言 ,x =(N-D) Fru, 这 导致 了 在 使 用 不 同 软 件 进行 分 析 时 ,研究 者 可 能 需要 针 
对 软件 的 计算 方法 对 公式 进行 微调 。 本 部 分 主要 展示 的 是 第 二 部 分 的 内 容 ， 即 已 知 So, df 和 a, 3 和 1-B 
已 知 其 一 时 如 何 使 用 G*Power 求 出 另 一 个 的 对 应 值 ,参考 前 述 对 G*Power 中 的 Generic y’ test 的 介绍 ， 这 
里 主要 分 三 种 情况 给 出 具体 的 操作 方法 与 示例 : 

(1) so=0,5 >0 

DE Al 1-B, 未知 61( 对 应 先 验 检验 力 分 析 ): 选择 检验 族 (Test Family) x? tests， 检 验方 式 (Statistical 
test) Generic y’ test, 检验 力 分 析 类 型 (Type of power analysis) 为 Sensitivity, 设 定 df. a 和 1-B, 求解 出 
51， 而 后 根据 Satorra—Saris 法 得 到 的 Fy 值 或 MacCallum 法 对 应 的 拟 合 指数 结合 补充 材料 C 部 分 的 公式 
推出 结果 。 

例 1: EA df= 51, 控制 a 水 平 为 0.05, 求 使 得 当 RMSEA 的 真实 值 不 小 于 0.08 时 正确 拒绝 零 假 设 
的 概率 不 低 于 95% 的 样本 量 No jan 

解 : 打开 G*Power， 选 择 检验 族 为 x tests, 检验 方式 为 Generic x? test, 检验 力 分 析 类 型 为 Sensitivity， 
设 定 4f= 51、a = 0.05、1-B = 0.95, 求 得 非 中 心 参数 为 43.400, 根据 公式 (C.2a), 求 出 样本 量 为 133.965， 
向 上 取 整 得 Norn = 134 Ao 

OBA ön AA 1-B( 对 应 后 验 检验 力 分 析 、 使 用 线性 插值 估计 样本 量 时 的 先 验 检验 力 分 析 ): 选择 
检验 族 为 x? tests, 检验 方式 为 Generic y test, 检验 力 分 析 类 型 为 Post hoc, E dfa 和 51, 求解 出 1-B。 

例 2: 同 例 1, 使 用 设 定 多 个 N 的 方式 结合 线性 插值 法 估计 Np。 

解 : 打开 G*Power， 选 择 检验 族 为 x? tests, 检验 方式 为 Generic x test, 检验 力 分 析 类 型 为 Post hoc, 
WE df=51、a= 0.05、 根 据 公式 (C.2a) 计 算出 N 为 S1, 101, 151 时 对 应 的 非 中 心 参数 分 别 为 16.32, 32.64, 
48.96, 分 别 代 入 , 求 得 对 应 的 1-B 为 0.566, 0.837, 0.975。 使 用 线性 插值 法 求 出 Npa AN 142 Ao 

(2) 80 = k(N-1) >0,61=0 

DEA 1-B， 未知 50( 对 应 先 验 检验 力 分 析 ): 根据 正文 5.6 部 分 的 结论 , 构建 与 其 对 偶 情 形 一 一 5。= 
0, 51 = KN 一 7) > 0， 而 后 交换 o 5 p, 选择 检验 族 为 x2 tests, 检验 方式 为 Generic y’ test, 检验 力 分 析 类 型 
为 Sensitivity, E df, a 和 1-B, 求解 出 51/， 进 而 补充 材料 C 部 分 的 公式 求 得 No 

例 3: 已 知 df=51, 控制 oa 水平 为 0.05， 进行 等 效 性 检验 , WE RMSEA, = 0.08, 求 使 得 RMSEA 的 
真实 值 为 0 时 正确 拒绝 零 假设 的 概率 不 低 于 80% 的 样本 量 Neptuno 

解 : 构建 对 偶 情形 : 5。= 0, 根据 公式 (C.2a) 有 51 = 0.3264* (N-1), 此 时 有 o = 0.2, 1-B’ = 0.95. FT 
FF G*Power， 选 择 检验 族 为 x tests, 检验 方式 为 Generic y test, 检验 力 分 析 类 型 为 Sensitivity, W df= 
51、Q=0.2、1-B=0.95, 求 得 非 中 心 参数 为 31.564, 进一步 求 得 Nou, AN 98 人 。 

QAM So, 未 知 1-B( 对 应 后 验 检 验 力 分 析 、 使 用 线性 插值 估计 样本 量 时 的 先 验 检验 力 分 析 ): 根据 
正文 5.6 部 分 的 结论 , 构建 与 其 对 偶 情 形 一 一 8, = 0, 8; =K(N-1)>0, 而 后 交换 a 5 p, 选择 检验 族 为 x 
tests, 检验 方式 为 Generic y test, 检验 力 分 析 类 型 为 Criterion, 设 定 df、a 和 61, 求解 出 1-B。 

例 4: 同 例 2, 使 用 设 定 多 个 N 的 方式 结合 线性 插值 法 估计 Npa 

解 : 构建 对 偶 情 形 : 59 = 0, 根据 公式 (C.2a) 有 61 = 0.3264 - (N-1), 此 时 有 1-8 = 0.95。 设 定 入 分 别 
为 51, 76, 101, 求 得 对 应 的 61 为 16.32, 24.48 和 32.64。 打 开 G*Power， 选 择 检验 族 为 x? tests, 检验 方式 
H Generic Xx? test, 检验 力 分 析 类 型 为 Criterion, 设 定 4f= 51、1-B?= 0.95, 将 8 的 三 个 值 分 别 代 入 , RK 
得 对 应 的 a Fy 0.614, 0.370, 0.180。 使 用 线性 插值 法 求 出 Nerm 约 为 99 人。 

(3) 89 = k(N-1) > 0, 8; = q(N—-1) > 0, 8) > 8, 

此 种 情形 下 问题 过 于 复杂 , 在 此 仅 给 出 一 种 情况 下 的 操作 : 

已 知 Sn 未 知 1-B( 对 应 后 验 检验 力 分 析 、 使 用 线性 插值 估计 样本 量 时 的 先 验 检验 力 分 析 ): 注意 到 ， 
给 定 自 由 度 时 中 心 化 x 分 布 的 右 尾 概率 等 价 于 其 临界 值 ， 因 此 可 以 以 此 为 “ 桥 ” 连 接 起 两 个 非 中 心 入 
分 布 。 首 先 选择 检验 族 为 x? tests, 检验 方式 为 Generic y? test, 检验 力 分 析 类 型 为 Criterion, 设 定 其 统计 
检验 力 为 等 效 性 检验 的 1-a, df, So, 求 出 中 心 化 好 分 布 的 右 尾 概率 ， 记 为 Y。 然 后 选择 检验 族 为 x? tests, 
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检验 方式 为 Generic y test， 检 验 力 分 析 类 型 为 Post hoc, 代入 a=y, 非 中 心 参数 为 61 和 自由 度 , 求 得 对 
应 的 统计 检验 力 (此 时 的 统计 检验 力 为 B 而 非 1-B)。 特 别 需 要 注意 的 是 ， 这 种 以 中 心 化 邓 分 布 为 “ 桥 ” 
的 方法 受 限 于 G*Power 对 a 的 计算 限制 , G*Power 允许 的 a 范围 是 1e-6 到 1-1e-6, 这 也 意味 着 此 种 情形 
下 G*Power 的 作用 是 有 限 的 ,对 于 此 种 情形 ,若非 中 心 ?分 布 的 自由 度 足 够 大 时 ,也 可 借用 正 态 分布 对 
非 中 心 ?分 布 进行 近似 从 而 求 出 统计 检验 力 的 近似 值 。 

例 5: 已 知 4f=51, 样本 量 为 301, WE RMSEA, = 0.08 进行 等 效 性 检验 , 结果 显示 RMSEA] = 0.043, 
p < 0.05, 求 统计 检验 力 。 

解 : 根据 公式 (C.2a) 计 算出 ôo = 97.92, 61 = 28.290。 打 开 G*Power， 选 择 检验 族 为 x tests, 检验 方式 
H Generic x? test, 检验 力 分 析 类 型 为 Criterion, 设 定 d= 51、1-B= 1-0.05 = 0.95, 入 = 97.92, R4 a X 
1.014e-6。 然 后 选择 检验 力 分 析 类 型 为 Post hoc, 设 定 4f= 51、 和 = 28.290, a = 1.014e-6, 求 得 统计 检验 
力 为 0.985。 

例 6: 已 知 df= 51, 样本 量 为 439, WE RMSEA = 0.08 进行 等 效 性 检验 , 结果 显示 RMSEA1 = 0.067, 
p <0.05, 求 统计 检验 力 。 

解 : 根据 公式 (C.2a) 计 算出 ôo = 142.9632, 5, = 100.2753。 打开 G*Power， 选 择 检验 族 为 x? tests, 检验 
FN Generic X2 test, 检验 力 分 析 类 型 为 Criterion, W df=51, 1-ß= 1-0.05 = 0.95, 4 = 142.9632, 求 
得 a 为 4.064e-012。 注意 到 解 得 的 a 远 超出 了 G*Power 允许 的 范围 , 但 同时 注意 到 本 例 中 df 较 大 ， 可 以 
考虑 使 用 正 态 分 布 对 非 中 心 ?分 布 进行 近似 ， 此 时 读 出 a 为 4.064e-012 时 对 应 的 中 心 化 x? Artn ie 
为 152.944， 同 时 计算 非 中 心 参数 为 8 、 自 由 度 为 51 的 非 中 心 x 分 布 的 均值 和 标准 差 .根据 非 中 心 x 
分 布 的 相关 结论 有 : 
MH=6+df =100.2753 + 51=151.275;6 =V24df +48 =V2x51+4x100.2753 = 22.430 (D.1) 

以 正 态 分 布 对 非 中 心 x 分 布 进行 近似 , 可 求 得 均值 为 151.275, 标准 差 为 22.430 的 正 态 分 布 割 点 为 
152.944 时 的 左 尾 概率 约 为 0.530， rae B 部 分 得 到 的 理论 值 (0.546) 还 是 较为 相近 的 。 


补充 材料 EE: 案例 中 Monte Carlo 模拟 部 分 各 参数 估计 情况 、 相 对 误差 以 及 95% 置 信 区 间 覆 盖 率 


REL 模型 拟 合 检验 先 验 检验 力 的 Monte Carlo 分 析 


模拟 样本 量 n= 101 模拟 样本 量 n = 201 
路 径 HEE 估计 ”均值 ”标准 gyp SE ”人 计 均值 ”标准 yeg 均值 标准 
均值 ”标准 误 误 均值 误差 RR ”均值 ”标准 误 误 均值 误差 ” 误 误差 


POS 一 Xl 0.870 0.867 0.083 0.084 0.953 0.004 0.018 0.870 0.059 0.058 0.952 0.000 0.009 
POS 一 X2 0.900 0.896 0.082 0.082 0.946 0.005 0.011 0.901 0.058 0.057 0.941 0.001 0.009 
POS > X3 0.750 0.745 0.088 0.089 0.947 0.007 0.008 0.752 0.063 0.062 0.952 0.003 0.013 
POS > X4 0.550 0.545 0.096 0.096 0.944 0.009 0.006 0.554 0.068 0.068 0.953 0.007 0.002 
DES > X5 0.590 0.589 0.099 0.095 0.956 0.001 0.035 0.590 0.070 0.069 0.958 0.001 0.013 
DES 一 X6 0.720 0.719 0.093 0.096 0.939 0.002 0.030 0.725 0.066 0.068 0.942 0.006 0.029 
DES > X7 0.580 0.577 0.091 0.096 0.943 0.005 0.049 0.581 0.065 0.066 0.943 0.002 0.017 
DES 一 X8 0.710 0.702 0.094 0.092 0.949 0.011 0.020 0.711 0.067 0.066 0.953 0.001 0.003 
DES > X9 0.590 0.587 0.091 0.092 0.947 0.006 0.019 0.593 0.064 0.066 0.946 0.005 0.025 
ANG 一 X10 0.730 0.731 0.093 0.093 0.946 0.001 0.006 0.733 0.066 0.066 0.947 0.005 0.010 
ANG — X11 0.780 0.779 0.093 0.094 0.946 0.001 0.013 0.780 0.066 0.069 0.934 0.000 0.040 
ANG > X12 0.570 0.574 0.092 0.095 0.937 0.007 0.032 0.572 0.066 0.066 0.948 0.003 0.012 
Xl XI 0.240 0.237 0.057 0.059 0.932 0.013 0.041 0.240 0.040 0.040 0.955 0.001 0.018 
X2 e X2 0.190 0.189 0.055 0.057 0.948 0.006 0.029 0.188 0.039 0.039 0.951 0.011 0.003 
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模拟 样本 量 n = 101 模拟 样本 量 n = 201 
BE HEA 估计 ”均值 ”标准 ggg BEO PE Mk H E ge 均值 标准 

均值 ”标准 误 误 均 值 OEO 误差 误 误差 ”均值 标准 误 误 均值 O 误差 误 误 差 
X3 过 X3 0.440 0.437 0.072 0.073 0.932 0.007 0.015 0.438 0.051 0.052 0.944 0.005 0.013 
X4 > X4 0.700 0.699 0.104 0.104 0.924 0.002 0.001 0.696 0.073 0.074 0.940 0.006 0.015 
X5 > X5 0.650 0.642 0.102 0.106 0.934 0.012 0.032 0.646 0.073 0.075 0.939 0.006 0.038 
X6 = X6 0.470 0.466 0.085 0.088 0.947 0.009 0.031 0.461 0.060 0.062 0.924 0.019 0.023 
X7 e X7 0.540 0.532 0.086 0.088 0.929 0.014 0.017 0.536 0.062 0.063 0.941 0.007 0.029 
X8 > X8 0.490 0.488 0.087 0.088 0.936 0.003 0.013 0.488 0.062 0.061 0.946 0.004 0.014 
X9 <> XO 0.520 0.516 0.085 0.086 0.937 0.007 0.014 0.519 0.060 0.064 0.933 0.002 0.057 


X10 e X10 0.450 0.441 0.086 0.091 0.921 0.021 0.055 0.445 0.062 0.063 0.940 0.011 0.022 
Xll e X11 0.390 0.384 0.085 0.087 0.950 0.015 0.018 0.391 0.061 0.062 0.949 0.003 0.012 
X12 e X12 0.550 0.539 0.088 0.091 0.929 0.021 0.031 0.546 0.063 0.062 0.939 0.007 0.014 
POS © DES 0.230 0.227 0.111 0.117 0.935 0.014 0.050 0.230 0.079 0.080 0.937 0.002 0.009 
POS ANG 0.400 0.398 0.103 0.110 0.929 0.004 0.067 0.401 0.073 0.074 0.953 0.002 0.004 
DES e ANG 0.780 0.779 0.071 0.072 0.938 0.001 0.016 0.781 0.050 0.050 0.943 0.002 0.000 


RE-2 模型 参数 检验 先 验 检验 力 的 Monte Carlo 分 析 


模拟 样本 量 n= 301 模拟 样本 量 n = 401 
路 径 。 理论 值 估计 均值 标准 gep 均值 标准 ”估计 均值 RE, ”均值 标准 
均值 “标准 误 误 均 值 “”” 误差 误 误差 ”均值 标准 误 误 均值 “” ” 误差 误 误差 


RES—SRG 0.180 0.179 0.057 0.057 0.951 0.006 0.001 0.181 0.049 0.048 0.949 0.001 0.026 
PSS—SRG 0.180 0.183 0.057 0.059 0.942 0.017 0.034 0.180 0.049 0.050 0.948 0.002 0.016 
ECI>SRG 0.180 0.180 0.057 0.056 0.948 0.002 0.016 0.179 0.049 0.050 0.942 0.006 0.013 
STR=SRG 0.038 0.038 0.053 0.053 0.949 0.010 0.009 0.040 0.046 0.046 0.950 0.036 0.000 
Ser>SRG 0.038 0.040 0.053 0.053 0.955 0.057 0.001 0.037 0.046 0.046 0.959 0.039 0.012 
Age>SRG 0.038 0.039 0.053 0.054 0.949 0.023 0.019 0.039 0.046 0.046 0.953 0.023 0.008 
PSS 一 RES 0.220 0.220 0.056 0.056 0.955 0.002 0.011 0.221 0.049 0.047 0.959 0.006 0.027 
ECI>RES 0.220 0.222 0.056 0.056 0.949 0.008 0.005 0.219 0.049 0.046 0.960 0.005 0.058 
STR=RES 0.047 0.047 0.054 0.057 0.940 0.002 0.049 0.045 0.047 0.048 0.947 0.031 0.018 
Ser 一 RES 0.047 0.046 0.054 0.054 0.949 0.010 0.006 0.046 0.047 0.048 0.947 0.010 0.017 
Age 一 RES 0.047 0.045 0.054 0.054 0.955 0.046 0.006 0.048 0.047 0.046 0.955 0.031 0.022 
ECI 一 PSS 0.282 0.282 0.055 0.054 0.961 0.001 0.017 0.283 0.048 0.046 0.958 0.002 0.036 
STR 一 PSS 0.060 0.060 0.055 0.056 0.949 0.010 0.008 0.062 0.048 0.049 0.951 0.038 0.020 
Ser—PSS 0.060 0.060 0.055 0.056 0.949 0.005 0.003 0.060 0.048 0.048 0.943 0.005 0.010 
Age 一 PSS 0.060 0.055 0.055 0.059 0.935 0.077 0.062 0.058 0.048 0.048 0.942 0.025 0.003 
n 为 501 和 601 时 的 参数 估计 情况 见 下 页 续 表 ， 加 粗 的 表示 该 路 径 系数 是 关键 参数 ; RES: 情绪 调节 自我 效能 感 , PSS: 
领悟 社会 支持 , ECI: 情绪 创造 力 , SRG: 压力 后 成 长 , STR: 疫情 相关 压力 , Ser: 疫情 严重 程度 , Age: 年 龄 
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R E-2( 续 ) ”模型 参数 检验 先 验 检验 力 的 Monte Carlo 分 析 ( 续 表 ) 


模拟 样本 量 n = 501 模拟 样本 量 n = 601 

路 径 EDE fhi 均值 ”标准 ey BO 标准 估计 均值 估计 均值 ”均值 标准 
均值 ”标准 误 误 均值 “” RA 误 误 差 均值 ”标准 误 、 均 值 ”标准 误 、 误 差 误 误差 

RES—>SRG 0.180 0.179 0.044 0.043 0.952 0.008 0.013 0.178 0.040 0.041 0.940 0.011 0.022 
PSS~SRG 0.180 0.180 0.044 0.046 0.943 0.003 0.035 0.182 0.040 0.041 0.944 0.009 0.023 
ECITSRG 0.180 0.183 0.044 0.043 0.952 0.015 0.019 0.180 0.040 0.040 0.941 0.001 0.000 
STRSRG 0.038 0.038 0.041 0.040 0.956 0.017 0.019 0.039 0.038 0.038 0.951 0.009 0.006 
Ser >SRG 0.038 0.039 0.041 0.041 0.954 0.028 0.002 0.038 0.038 0.039 0.947 0.005 0.031 
Age>SRG 0.038 0.041 0.041 0.041 0.957 0.064 0.013 0.039 0.038 0.037 0.947 0.015 0.009 
PSS 一 RES 0.220 0.221 0.044 0.042 0.962 0.003 0.037 0.218 0.040 0.041 0.951 0.010 0.018 
ECI 一 RES 0.220 0.221 0.044 0.042 0.958 0.006 0.045 0.220 0.040 0.039 0.952 0.002 0.029 
STR 一 RES 0.047 0.050 0.042 0.043 0.947 0.067 0.020 0.047 0.038 0.038 0.950 0.005 0.019 
Ser—>RES 0.047 0.045 0.042 0.041 0.958 0.029 0.010 0.046 0.038 0.039 0.943 0.017 0.009 
Age—RES 0.047 0.047 0.042 0.042 0.947 0.013 0.003 0.048 0.038 0.038 0.955 0.038 0.005 
ECI—>PSS 0.282 0.282 0.043 0.042 0.948 0.001 0.010 0.280 0.039 0.037 0.963 0.008 0.049 
STR 一 PSS 0.060 0.059 0.043 0.042 0.956 0.018 0.023 0.061 0.039 0.040 0.947 0.012 0.011 
Ser 一 PSS 0.060 0.059 0.043 0.043 0.946 0.010 0.010 0.061 0.039 0.041 0.945 0.012 0.041 
Age 一 PSS 0.060 0.060 0.043 0.043 0.945 0.000 0.000 0.059 0.039 0.040 0.948 0.009 0.029 


